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We present results for the static inter-quark potential, lightest glueballs, light hadron 
spectrum and topological susceptibility using a non-perturbatively improved action on a 
16 X 32 lattice at a set of values of the bare gauge coupling and bare dynamical quark 
mass chosen to keep the lattice size fixed in physical units (~ 1.7 fm). By comparing 
these measurements with a matched quenched ensemble, we study the effects due to two 
degenerate flavours of dynamical quarks. With the greater control over residual lattice 
spacing effects which these methods afford, we find some evidence of charge screening and 
some minor effects on the light hadron spectrum over the range of quark masses studied 
{Mps/Mv > 0.58). More substantial differences between quenched and unquenched 
simulations are observed in measurements of topological quantities. 



I. INTRODUCTION 

Over recent years, considerable effort has gone into probing QCD beyond the quenched approximation. 
For recent reviews see [0-0 and for results using a different improvement scheme see B. Because of 
the impressive agreement of the quenched approximation (see e.g. B) with experiment for the spectrum 
and other easily accessible quantities, the effects of dynamical quarks in these are expected to be quite 
small. It is difficult to isolate physical effects which are unambiguously due to their inclusion, in part 
because of the need for high statistics. On currently available machines this requires coarse lattices. 
The use of 0{a) non-perturbatively improved fermions has been suggested as a means of controlling and 
reducing discretisation errors [Q. In an earher paper Q, first results of the UKQCD Collaboration using 
a preliminary value of the improvement coefficient Csw were presented. It was found that the effective 
lattice spacing, as measured by Sommer's intermediate scale parameter tq pj, depended quite strongly on 
the bare quark mass at fixed gauge coupling. However, the effect of dynamical quarks on easily accessible 
physical observables was very weak and difficult to disentangle from those induced by other changes in 
the simulation parameters. Eventually, one might hope to perform detailed studies over the full space 
of parameters including bare gauge coupling, quark mass(es) and lattice volume. In the meantime, less 
ambitious studies may still serve as a guide to those regions of parameter space where physical effects 
may be found. 

In this paper, we present results of further simulations over a range of sea quark masses. For these 
simulations, we have used the final published values of the 0{a) improvement coefficient Csw Q] ^^id have 
attempted to reduce variations due to residual discretisation errors and finite volume effects by working 
at fixed lattice spacing. In order to achieve the latter, we have used matching techniques described in 
an earlier work llfl] to help obtain ensembles of configurations whose lattice spacings, as defined by the 
scale To, are as closely matched as practicable. We present results for the spectrum and potential on, 
or close to, a single fixed rg trajectory in the (/3, k) plane which extends from quenched configurations 
(k = 0) to the lightest accessible sea quark mass. We choose tq to set the scale since it has no valence 
quark complications and is determined by intermediate scale properties of the static potential. These 
properties are expected to be less sensitive to charge screening (short range) and string-breaking (long 
range) effects arising from dynamical light quarks. 

We interpret our results in the spirit of partial quenching. That is, we study chiral extrapolation in the 
valence quark masses of light hadron masses using both quenched and partially unquenched configura- 
tions. We find that, with the available statistics, the quality of these valence extrapolations is uniformly 
good. By studying the spectra so obtained we search for evidence of the influence of light dynamical 
quarks. We also study the behaviour of the topological susceptibility in the presence of dynamical quarks. 
Our data sample includes measurements made with equal valence and sea quark masses. 

The plan of the rest of the paper is as follows. Section || contains brief details of the simulation methods 



and parameters. In section III, we review the matching techniques used to set up simulations at similar 
lattice spacings. We present results in section |^ for the static potential in QCD and use it to define a 
lattice scale. In section M, we present results for the light hadron spectrum including some measurements 
of the lightest glueball masses. Section VI contains results from topological charge and susceptibility 



measurements. Finally, our conclusions are summarised in section VII 



Some preliminary results from these analyses have been presented elsewhere 11 16 



II. SIMULATIONS WITH IMPROVED WILSON FERMIONS 



Details of our implementation of the Hybrid Monte Carlo simulation algorithm |0] and its performance 
can be found in our earlier paper S. Here, we summarise for convenience some key features. For the 



lattice action we used a standard Wilson action for the gauge fields together with the Sheikholeslami- 
Wohlert 0(a)-iniprovcd Wilson gauge-fermion action ||l8|] : 



where 



5[C/,V,V'] = 5G[C/] + 5f[f/>,^], (1) 



Sg [U] =pWa^PY.(^- ^R° Tr Up) (2) 



and 

Here, Up is the usual directed product of gauge link variables and S^ is the standard Wilson fermion 
action, 

Sf =Y.^{x)ij{x) -kY, (?(a;)(r - ^^,)U^{x)^{x + fi)+ Jj{x + fi){r + j^)Ul{x)i:{x)) (4) 

with the Wilson parameter chosen as r = 1. The spin matrix is cr^^ = ^ [7^^,71/], and F^^{x) is the field 
strength tensor 

F^.{x) = i (/^.(x) - fl,{x)) (5) 

where /pi/(x) = U^i,{x) + U,j^^^{x) + {/_iy__^(a;) + t/_^^^(a;) is the sum of four similarly oriented (open) 
plaquettes around a site, x Q. 

Beyond tree level, the improvement coefficient Cgw is a function of the gauge coupling /3(= 6/5^). In 
the studies reported here, we have used those values determined non-perturbatively by the Alpha Col- 
laboration and summarised by an interpolation formula M. For example, at /? = 5.20 we have used 
c«w = 2.0171 0. 

We have used two degenerate flavours of dynamical quarks in these simulations. The bare quark mass is 
controlled by the hopping parameter k,. Restoration of (spontaneously broken) chiral symmetry, requires 
extrapolation in n to the critical value Kcrit at which the pion is effectively massless. As discussed above, 
we will often discuss the situation encountered in the quenched approximation, where the dynamical 
(sea) quark mass parameter {k^'^'^) is fixed (at in the quenched case) while the chiral extrapolation is 
performed in the valence mass parameter (k™') only. This is often referred to as a partially quenched 
approximation. It is particularly relevant where the dynamical quark mass is still quite heavy and where 
there is no realistic prospect of approaching the (degenerate) light quark chiral limit in both parameters. 



A. Simulation parameters 

Since these simulations were the first to be done on a reasonably large lattice (16'^ x 32) using the fully 
improved value of Csw, there was little guidance available on the choice of simulation parameters. We 



* Although the effect of 0{a) improvement is not expected to be as sensitive as the quoted number of significant 
figures suggests, the action and lattice observables do depend quite strongly on this parameter. For reasons 
of reproducibility we have therefore used a 4 decimal place representation of the Csw formula in generating 
configurations. 



chose /3 — 5.20 as the lowest value at which a reliable value of Cgw was available [Q. The aim was to 
obtain as large a physical volume as practicable with the available computing resource. The use of an 
improved action was expected to offset (at least partially) the relatively coarse lattice spacing which this 
implied. Equilibration was carried out through a sequence of dynamical quark masses: k'*"^ = 0.13000, 
0.13350, 0.13400, 0.13450 to 0.13500. The first production run was then carried out at k"'=^ = 0.13500 
starting at trajectory number 10010, where trajectories were of unit length. Configurations were stored 
after every 10 trajectories although a larger separation was used for most operator measurements (see 
below). 

Further simulations at higher quark masses (k*"^^ = 0.13450, and 0.13400) and slightly shifted (3 were 



then performed. The shifts in (3 were estimated using the methods described in section III and were 
designed to maintain a constant lattice spacing as defined by tq. 

To complete the comparison of unquenching effects, we performed pure gauge simulations using a standard 
update algorithm, heat-bath with over-relaxation. Again, the /3 value was chosen to keep tq at the value 
measured on the ensemble obtained at {P,k^°^) — (5.20,0.13500). The only exceptional configuration 
found was within the quenched configurations, and this was only apparent for one of the k™' studied 
This configuration was excluded from further analyses. 

An additional substantial, but unmatched, simulation was then performed at (/3, k^°^) — (5.20, 0.13550). 
This ensemble of configurations was analysed along with the matched ensembles providing further infor- 
mation on behaviour at light quark mass. A simulation at even lighter quark mass (k'^"^ = 0.13565) was 
begun. Where relevant, some preliminary results are presented here. Table | contains a summary of the 
run parameters for each ensemble. 

The bulk of the simulations were carried out in double precision. This followed initial concerns over the 
effect of rounding errors on reversibility. Detailed analysis of these and related effects have been carried 
out and have been reported elsewhere [|l9| . This work shows that, at least for present volumes and step 
lengths, the algorithm is reversible and stable for all practical purposes, even when implemented in single 
precision. 



B. Autocorrelations 

We made autocorrelation measurements from the average plaquette value measurements on every trajec- 
tory. The methods used were those described in detail in our earlier paper l^]. As shown in Table ||, 
the observable autocorrelation (from the plaquette) is of order 20 and so we have adopted a separation 
of 40 trajectories as standard in the analysis which follows. Nevertheless we keep in mind that subtle 
longer-term autocorrelations, not directly measurable, may still be present and so we have done additional 
checks on our statistical error estimates by re-binning the measurements. In the present data sample, we 
have not found any evidence of such correlations. 

Further measurements of the integrated autocorrelation time have been attempted for the potential 



(section |VD and the scalar glueball (section |V G| ) . At the lightest quark mass {k^'^'^ = 0.13565) autocorre- 
lations were estimated from effective mass (potential energy) measurements made every 20 trajectories at 
various lattice distances {r/a = 1 — 5) and Euclidean times {t/a = 3 — 5). The measured integrated auto- 
correlation times varied from 10 to 20 trajectories with large errors (typically ±8). For the scalar glueball, 
the integrated autocorrelation time for effective masses was in the range 25 to 30 at k^''^ = 0.13500 and 
0.13550. 

It is noteworthy that the autocorrelation is significantly less in the current simulations than in our 
previous runs |q| at comparable quark masses but different Csw The current simulations use the fully 
non-pcrturbatively improved value of Csw It is further noted that r'"* appears to decrease, if anything, 
with decreasing quark mass. This is contrary to the simple expectation that, as the lattice correlation 



length (typically given by the inverse pion mass) increases, then so should the correlation in computer 



time. A similar effect is evident in the decorrelation properties of the topological charge (see section |V|) • 
It is possible to reproduce such behaviour in simple models. The integrated autocorrelation time, which 
determines the size of the errors, can decrease even in the presence of increasingly long correlation modes 
simply due to increased noise induced by dominant short correlation modes. 

To illustrate this point, consider first the following simple model consisting of a single Markov chain 
x{t),t = 0,1, 2,... 

x{t) ^ ax{t - 1) + z{t) , a;(0) = 0, (6) 

where the z{t) are uncorrelated Gaussian noise of unit variance and < a < 1. It is simple to show that, 
for sufficiently long chains. 



p^{t) = a* = e-*/""'(") and so (7) 

1 l + g 
21-a 



r-P(x) = -lna, r'"* (x) = ^ i±^ . (8) 



Here Px(t) is the normalised autocorrelation function for the observable x. The corresponding results for 
finite length chains are also calculable so one can study the effects of using limited statistics to estimate 
autocorrelation times. Here we stick to the infinite chain approximation. For a ^ 1 (i.e. for large 
r-'^P(x)), 

r'"*(x) = r'='^P(a;) + l/(12r'='^P(a;)) + C'((t'=^p)-3) . (9) 

Real HMC data for p{t) do not, of course, show a simple exponential behaviour and so it is useful to 
consider the next simplest model which contains two independent correlation modes with relative coupling 
strength r: 

X{t)^xi{t)+rx2{t), (10) 

Xi{t) = a,Xi{t - 1) + Zi{t) (i = l,2). 

The integrated autocorrelation time for X{t) is given by 

r'"*(X) - 7yiT'"'(a;i) + ??2t'"*(x2) , where (11) 

Thus the relation between the integrated autocorrelation time and the actual correlations present in Xt 
is no longer straightforward. There may be quite long correlations present {t°^^{x2) ^ t°^p{xi)) but, 
depending on the relative strength of the modes (given by r) , the 'weighted average' represented by the 
above formula can give a result bearing no relation to either r°''P(xi) or t°^'^{x2). 

The possibility of such behaviour makes it essential to check decorrelation for individual observables 
explicitly using binning techniques. 



C. Finite-size effects 



In retrospect, the value of k^°^ = 0.13500 turned out to be somewhat conservative, in that the correspond- 
ing ratio of Mps/My is quite large (0.70, see section Q). The choice was based on preliminary estimates 
of the limiting algorithm performance and on measurements of the effective lattice spacing as described 



in section IV . It was felt that decreasing the quark mass further would decrease the effective lattice size 



to a point where finite-size efi^ects would become a problem. In our earlier analysis of finite-size effects 
(at least as far as they affect the potential and light hadron spectrum for Mps/Mv^O.7) we found that 
such effects were negligible provided 

L/ro>3.2. (13) 

This corresponds to a spatial extent of around 1.6 fm and is satisfied by all but our lightest quark mass 



data set, as shown in Table III. Further investigations may be called for, given the concerns over the 
baryon mass spectrum noted in [^ (see sec.VC). 

III. MATCHING SIMULATION PARAMETERS 

In a previous paper [Q we have described techniques which allow one to use unbiased stochastic estimates 
of the logarithm of the fermion determinant to determine, approximately, curves of constant observable 
in the space of simulation parameters. 

A. Determination of fixed observable curves 



The approximate character of the formalism arises from two sources. First, the log of the fermion 
determinant is only determined stochastically on each configuration and the corresponding fluctuations 
are proportional to the lattice volume. Second, a linear approximation is used when dealing with small 
changes so that these curves may only be determined locally. In the present application, the parameter 
space of interest is the (/3, k) plane and the observable of interest is the QCD static potential scale 



parameter tq (see section IV). 



To first order in small parameter changes (5/3, 5k), the shift in the lattice operator F is given by |l3| 

(SF) = [{FWn) + (^^>c.w]<5/3 + {F^)^^^- (14) 

The quantity 



Csw 



dp 

is well determined 0] and so the identification of constant F curves 

{5F) = (15) 

reduces to measuring correlations of the form 

{FWa) and {FSf) . (16) 



Here, (A) denotes the connected part (^4 — {A)) of the operator A. We refer readers to |10| for a 
detailed discussion of the stochastic evaluation of T = TrLnM^M. Here, M is the fermion matrix 
including the 0{a) improvement term. The methods are based on a Lanczos implementation of Gaussian 
quadrature [E0| . Recent progress in understanding the nature of roundoff errors in the finite arithmetic 
Lanczos process assures us that this application of the Lanczos process, unlike the basic algorithm itself, 
is highly stable with respect to roundoff [pl| . 



B. Matching tq 

Detailed tests of the matching procedures have been carried out using the average plaquette, which is 
very accurately measured, and a variety of Wilson loops llG]. Some tests using ro were also carried out 
successfully on modest-sized lattices. The present work represents the first application, in earnest, to 
production-size lattices. Since the fluctuations in T are extensive quantities, we expect there will be a 
limit on the size of lattices where usefully accurate matching estimates may be made with a given amount 
of work. The work required has been analysed in some detail in |lO| . 



The correlations ( [16D require measurements of F on each configuration. These are available for operators 
such as Wilson loops but not for physical quantities such as hadron masses and rg . Rather than determine 
the fixed rg curve directly using dlq), we use (M) to estimate the required gauge correlators at nearby 
points in parameter space. We then extract the potential, and hence rg/a, at the nearby parameter values 
from these 'shifted correlators'. This allows estimates of the partial derivatives with respect to (3 and k 
and hence the shift SP^ required to compensate a particular change in k^°^. 

dF dF 

where, in the present application, F — tq ^ ro/an 

Using an ensemble of 100 configurations at {(3,k) — (5.2,0.1350) for all correlator measurements, we 
estimated that a shift of 

(5/3 = 0.057 ± 0.033 

would be required so as to match the value of fg at (5.2 -I- Sp, 0.1345) with that at (/?, k) = (5.2, 0.1350). 
A simulation run at (/3, k) = (5.26, .1345) confirms that rg, and hence the effective lattice spacing, is 
indeed well matched (see Table 0) . 

However, it is clear from the size of the statistical errors that estimates of 6/3 obtained on these config- 
urations cannot be relied upon, in general, to predict matched parameters with great accuracy without 
further checks. The level of success in achieving fg matching can be gauged from Table M. 

The above methods for matching parameters are only applicable for small shifts. To obtain the shift for 
matching quenched simulations {Sk = —k) we have used tabulated values 122] of fg (lattice spacing a) to 
provide an initial estimate. Since there are systematic differences arising from slightly differing methods 



for extracting the potentials and fg (see section IV), we used this only as an initial guide. Following 



direct measurement of fg with our own techniques, we then made a further small shift in /3. The results 
are shown in the next section. 



IV. THE QCD STATIC POTENTIAL 

We have determined the static inter-quark potential V{r) using standard methods and used it to search 
for signs of charge screening and string breaking, as well as to determine the physical scale. 



^Here, and in what follows, we use the notation A to denote a physical quantity A expressed in lattice units. 



A. Extraction of the potential 

The methods foUow those originally proposed by Michael and collaborators p3,pj| . A variational basis of 
generalised Wilson loops is constructed from gauge links which are 'fuzzed' in the spatial directions 123]. 
The spatial paths between the static sources include a limited number of off-axis directions as well as 



those along the lattice axes (see the lower half of Table IV). A transfer matrix formalism is then used to 
extract the Euclidean time energy eigenstates which are related to solutions of the generalised eigenvalue 
equation 

M/„-(r,i)0(r)f) = AW(r;t,io)W^.j(r,io)0(r)f , t,j,k^O,l. (18) 

Here, we have used two levels of 'fuzzing' so giving a 2 x 2 eigenvalue equation. We used level (unfuzzed) 
and level 16 which means sixteen transformations of the spatial links. The link/staple weighting used 
was 2.5. This choice of fuzzing parameters was made so as to give a satisfactory variational basis with a 
modest amount of computational effort. Initial tuning experiments were made using 20 configurations at 
(/3, K^'^'^) = (5.2, 0.13500) and repeated on a corresponding matched ensemble of quenched configurations. 
Expanding the basis to three levels of fuzzing did not significantly improve the resulting effective mass 
plateaux extracted as described below. 

In principle, one could use the largest eigenvalue A'-''-'(r;i, tp) for large i,io to estimate the potential. In 
practice, however, the eigenvalue system becomes unstable at large t, particularly when modest numbers 
of configurations are used, as is often the case in dynamical fermion studies. Instead, we used the leading 
eigenvector (j){r)^"\ corresponding to \^°\r;t,to) at i = 1, to = 0, to project onto the approximate 
ground state |§,^. The resulting correlator Wo{t)was then used to form effective mass estimates for the 
approximate ground state 

LWo(r,t+l)J 

The ratio of the first two transfer matrix eigenvalues 

i?i,o = A(i)(r;l,0)/A(o)(r;l,0) (20) 

was used to help obtain improved estimates of the ground state energy with reduced contamination from 
the first excited state. To do this, the correlator Wo{t) was modelled as a sum of of two exponential 
terms 

W^o(r-,t)cx[AW(r;l,0)]* + ei,o[A«(r;l,0)]*. (21) 

One can easily show that, provided the contamination from the first excited state (ei,o) is small, the true 
ground state energy in such a model is given by 

Eoir) = -lnA(«)(i?,o) « ^oir,t) -R^r^t ^ 1) ^^^^ 

1 — -Ki,o 

Rather than search for plateaux in this quantity, we used a weighted mean of values from tmin to tmax 
where the weighting was inversely proportional to the statistical error (estimated via simple jackknife). 
To obtain the final quoted values we used (tmin, tmax) = (4, 5). In all cases, the difference due to increasing 
or decreasing the cut-offs by one was less than the statistical errors quoted. Overall statistical errors were 
estimated by bootstrap sampling. 

We also studied double exponentials fits to the effective mass Eo{r,t) using time-slices up to t = 8 and 
exponential fits to the full 2x2 matrix correlator. The fits, where stable, yielded results compatible with 
those obtained by the above methods. 

In Figure |l|, we show examples of the effective mass and corresponding extrapolated energy (|2^) used to 



determine V{r). The lattice potential values are collected in Table IV 



B. Determination of ro/a 

The potential V{r) can be used to determine the force between a static quark-antiquark pair separated 
by a distance r = |r| and hence to extract the Sommer scale parameter tq. This is a characteristic 
scale at which one may match the inter-quark force with phenomenological potential models describing 
quarkonia M . Specifically, it is defined by the solution of the relation 

rl^l =1.65. (23) 

Physically, vq ~ 0.49 fm. and we adopt this latter value when physical units are required. This definition 
of the physical scale has the advantage that one needs to know the potential only at intermediate distances. 
An extrapolation of the potential to large separation, which is conventionally performed to extract the 
string tension, is thus avoided. Hence, the procedure is well-suited to the case of full QCD for which the 
definition of a string tension, as the limiting value of the force, is not applicable. The string is, of course, 
expected to break at some characteristic distance rf,. 

Our determination of ro/a follows the procedures originally described in |2^ and recently adapted to 
provide a comprehensive study of the scale parameter in quenched QCD [E2|. That is, we perform fits to 
the parametrisation 

nr) = K,+ar-e[i]+/([i]-i) (24) 

where [j;] is the tree-level lattice Coulomb term 



■ 1, _ r (Pk cos(fc 

^rJ ^^i_ (2^ 4E-.iSin2(fc,/2) 



I- in / a K cosiK • r I , , 

[-]=4W 7^ .v^3 L2.. ,.. ■ (25) 



The parameter / is introduced so as to model further lattice corrections beyond tree level. We find 
that, for the coarse lattice spacings considered in this work, a tree-level parametrisation gives a poor 
description of the data for r ^ tq. 

Following [g2[, we use fits of the form(p4|) to provide a good description of the intermediate range potential. 
We then identify the fitted parameters as reliable estimates (up to 0{a?')) of the corresponding continuum 
version which, from the definition of tq, satisfies 

arl + e^c= 1.65 (26) 

and hence we extract our estimate of ro as 



ro = A/ . (27) 



Both on and off- axis measurements of the potential were used (see Table IV). We confirm the observa- 
tion ||22| that the value of ro extracted in this way is remarkably insensitive to changes in the fit range 
used. The individual parameters such as e and / are, however, quite sensitive. The point at r = (a, 0,0) 
was omitted from all fits since its inclusion was found to give an unacceptably high contribution to x^. 
The inclusion of data at the largest \r\ (> 8a) played little role in the determination of rg. Since a limited 
range of r is used to determine the parameters of (p4), one should treat the value of a with some caution. 
It does not represent a careful determination of the string tension which of couse is a large distance 
property and, strictly speaking, only meaningful in the heavy sea-quark limit. 

We present a summary of the results for ro in Table [V|. The systematic error estimates (shown as (+x— j/)) 
were determined by variations in the fitting range used for r and in the number of parameters used in 



the fit. The central values quoted were obtained using all potential data satisfying y^ < r < 8. As 
described in the next subsection, a term proportional to 1/r^ was tried. The systematic error estimates 
also include the effects of varying t^^i^ by one unit in the evaluation of the potential (see above). It is 
seen that, for the ensembles at (5.20,0.1350), (5.26,0.1345) and (5.93,0), the matching in fo (and hence 
in effective lattice spacing) is very good (well within statistical errors) while that at (5.29, 0.1340) is only 
slightly off (just over one standard deviation). The unmatched simulation at the lightest quark mass has 
a significantly smaller lattice spacing (seven standard deviations). 



It is worth noting (Table IV) that the absolute values of the potential are not matched even when tq is. 
The same is true for the average plaquette and the generalised Wilson loops themselves which go into 
the potential determination. All of these loop operators have large ultraviolet-sensitive contributions. 
In section M, we will comment further on the extent to which matching is observed in other physical 
quantities. 

The value fo = 4.714(13) for the quenched measurements at /3 = 5.93 may be compared to previous 
high statistics measurements in quenched simulations. The interpolating parametrisations of |22f| and ||28|] 
respectively suggest 4.757 and 4.741(18) in fair agreement with, but slightly larger than, our determination 
of this quantity at this particular value of f3. The slight discrepancy of our result with that of rcf. [27] 
amounts to about one standard deviation.' 

The JLQCD collaboration have presented preliminary results from an Nj = 2 simulation using the same 
action as the present work at (3 — 5.2, Csw — 2.02 and n = 0.1350,0.1355 [2^. The values of fg in this 
case are slightly smaller than those presented in Table 0. Note that the value of Csw used by JLQCD is 
very slightly different from ours. The methods used to extract the potential and fg apparently have much 
in common with those described above but we have not been able to check all the details. In particular, 
the errors so far presented by JLQCD are statistical only. 

As mentioned above we have used both on-axis and off-axis Wilson loops in our determination of rg. How- 
ever, different spatial orientations of Wilson loops differ by lattice artefacts of order a^. Thus, if on-axis 
loops are used exclusively to extract rg, then the result may not be consistent with a determination using 
other orientations, provided that the statistical accuracy is large enough to expose these discrepancies. 
For our Nf — 2 simulations the level of precision is about 1%, so that any significant discrepancy in 7'o 
due to different orientations will be hard to detect. In future high-statistics simulations with dynamical 
quarks, a cleaner procedure might be to define rg consistently for one particular orientation and to ex- 
tract rg from local interpolations of the force between static quarks. This is the approach used in refs. 
[9,25,27,8]. It has also been used for some of the ensembles presented here and for rg itself makes little 
difference (within the statistical errors). 



C. Charge screening 

In Figure we plot the static potential in units of rg. The zero of the potential has been set at r = rg. 
Overall, the presence of dynamical fermions makes little difference when plotted in physical units. The 
data are apparently well described by the universal bosonic string model potential Q which predicts 

[V{r) - V{ro)]ro - (1.65 - e) ('i- - 1 V e (^ - l) . (28) 

Of course, the fact that the scaled potential measurements all have the same value and slope at r = rg 
simply reflects the definition of rg . In Figure ^ we show the deviations from this model potential. Here 
e — tt/12 [|3y]. We note the following points: 

• At the shortest distances (see the points where |r| < 0.5rg) there are indeed deviations from the 
string model. 



The large fluctuations as a function of r jr^ indicate strong violations of rotational symmetry (see 



Table |V| for a list of separations used). 



• 



There is some slight evidence that the deviations depend systematically on the quark mass - compare 
the quenched points with those for the lightest values of k^ 



-sca 



• The matching of the data ensembles allows a clean comparison of the data at different quark masses. 

• There are no indications of string breaking, but we note that the distance probed at light quarks 
masses is not large, (r < 1.3 fm at the lightest quark masses used.) 



As discussed above, the parametrisation (p4) is not particularly efficient at describing the short-range 
interactions on the lattice. This is the case even though it allows, in a model-dependent way (/ ^ 0), 
for lattice artifacts beyond those expected at tree level. The fits for the effective charge e and associated 
parameter / are therefore sensitive to the fit range and any variation in the parametrisation. For example, 
we also considered allowing a term proportional to |r|~^ in an attempt to describe better the short distance 
potential. However, the coarseness of the lattice and crudeness of the parametrisation prevented reliable 
fits. In the continuum limit one would expect the short range potential to behave as 

where \x is some scale. Lowest-order perturbation theory then suggests an enhancement of some 14% in 
as arising from the change in the factor 33 — 2A'^/ when unquenching the theory (at fixed scale). Using 



the above parametrisation (24) we can see if such an effect is reflected in a corresponding increase of 18% 
in the parameter e. 

We have performed correlated fits to the potential with a constant choice of parametrisation and fit 
ranges. Some reasonable variation in the latter was then used to give an estimate of systematic errors. 
The fits for the central values of parameters included all data from Table IV satisfying \f2 < \r\/a < 9. 



The statistical errors were produced via an overall bootstrap of the full analysis (with 500 bootstrap 
samples). The results are included in Table M The coupling parameter e does seem to show an increase 

due to unquenching. For the matched ensembles the increase is 18 _]^o% in going from quenched to 
K^t^a ^ ,13500. 

Similar findings in the case of two fiavours of Wilson fermions have been reported by the SESAM-T^L 
collaboration |Q where an increase of 16 — 33% was found. 

For comparison with other scale determinations, we have included the fit parameter ^/a expressed in 
units of MeV as deduced from tq = 0.49 fm. We repeat the caveat offered above that the parameter ^/a■ 
reflects the medium-range shape of the potential and does not represent a definitive determination of the 
asymptotic string tension. Phenomenological models for the hadronic string suggest a value of around 
440 MeV. The energy scale determination based on ro/a is therefore some 6 to 7% higher than that based 
on the string tension. In the next section we compare the above scale determination with values deduced 
from the vector meson mass. 

Recently the MILC collaboration ||3^,Q has presented results of a comparison of the quenched static 
potential with that due to three flavours of staggered fermions. As in the present analysis, the authors 
have noted the strong influence of the dynamical quarks on the effective lattice spacing and have compared 
the shapes of the potential measured on matched ensembles. 



V. LIGHT HADRON SPECTRUM 

Throughout this section one of our main aims will be to uncover any unquenching effects in the light 
hadron spectrum. Because we have a matched data set, any differences can more directly be attributed to 



unquenching effects. However, the task of identifying differences is likely to be hard for those quantities 
which are primarily sensitive to physics at the same scale as that used to define the matching trajectory 
in the {/3, k^'^'^) parameter space (rg in this case). This is expected to be the case for the hadron spectrum 
considered here where the quark masses are still relatively heavy. 

Two-point hadronic correlation functions were produced for each of the datasets appearing in Table |. 
The interpolating operators for pseudoscalar, vector, nucleon and delta channels were those described 
in [Q. Mesonic correlators were constructed using both degenerate and non-degenerate valence quarks, 
whereas only degenerate valence quarks were used for the baryonic correlators. 



The hadronic masses are presented in Tables |VI| - |IX| . These are expressed in both lattice units (M = Ma), 
and in the dimensionless form roM. Note that the errors displayed are statistical only. We estimate that 
the systematic errors arising from different choices of fitting procedure are similar in size to the statistical 
errors. 

In the following, we review the main fitting procedures which were used to obtain the light hadron 
spectrum results. Further details of the fitting procedure can be found in |M. 



A. Fitting procedure 

We used the fuzzing procedures of |^ to generate correlators of the type LL, FL and FF where F denotes 
fuzzed, and L local operators. Conforming to our usual convention, FL means fuzzed at the source and 
local at the sink. The fuzzing radius was set to Rfuzz — 2. 

Effective mass plots for the three types of fuzzed correlators (LL,FL and FF) are shown in Fig. ^ for the 
f3 — 5.2, K^°^ — K™' — 0.13500 data set. Note that all the effective mass plots approach their asymptote 
from above. The FF correlator exhibits the fastest approach. This behaviour is universal throughout all 
the datasets. For technical reasons, the fuzzing procedures used in practice for the hadron correlators 
introduced some unbiassed stochastic noise. We have checked that this has indeed had no significant 
effect on the hadronic quantities presented here but has resulted in increased error estimates at the level 
of less than 10% for the pion and less than 20% for the nucleon. 

Correlated fits were used throughout the fitting analysis of the correlation functions and the eigenvalue 
smoothing technique of ||3^ was employed. Ensembles of 500 bootstrap samples were used to estimate 
the errors [p8[ . 

We performed a factorising fit which we now describe for the baryonic case. The three fuzzed correlators 
LL, FL and FF are fitted together, where the fitting function used for, say, the FL channel is 

Zo^Ztfe^""* + Zf Zfe""^*, 

and the LL and FF fitting functions are similarly defined (see e.g. pa). Note that both the coefficients, 
Zq^i and the masses mo,i are common to all the channels, and that the x^ comprises the individual x^ 
of the three channels and includes the correlation between different times and channels. 

For the mesonic case, we modify the above as usual by including the backward-propagating state, i.e. 
Q-mt _^ g-m(T-t) -^j^gj-g y jg ^j^g temporal extent of the lattice. 

Within these three different fitting types, a sliding window analysis was used to determine the optimal 
fitting range (tmin — ^max) [pOJ- In this analysis, fits for various imin were obtained with imax fixed 
generally to 15. Stability requirements in the baryonic sector forced imax = 14 in some cases. The masses 



so obtained are displayed in Tables VI -IX 



B. PCAC Mass 

The PCAC mass can be defined using tlie relation 

df_,Af,{x) = 2mpcAcP{x), 

where P{x) and A{x) are pseudoscalar and axial current densities. On the lattice, the following expression 
can be used to obtain an estimate of mpcAC [0 ■ 

, BaCa^p^ (0, t) + acAdldiCpp^ (0, t) , 

"^'"'" ^ ^ ^^^7;^) ^ 

= (r(i))+CA(s(i)) (30) 

where d^ is the temporal lattice derivative averaged over the forward, 9, and backward, 9*, directions, 
and (. . .) represents averaging over times, t, where the asymptotic state dominates. The correlators C 
are defined in Q . The value of the coefficient used is 

CA = -0.007565^, (31) 

with gQ = 6//3 (the bare coupling). This is the one- loop, dynamical value [|2[, and hence eq.(|30|) suffers 
from 0{a) errors. Table M shows the results for rripcAC for all the datasets with ca defined as in eg. (|3l|) . 



There has been some recent debate in the literature regarding the most suitable non-perturbatively 
improved ca value (see e.g. ||43|-|45|]) and a reliable value may not yet have been determined. In the 
absence of a non-perturbatively improved value of ca (for Nj = 2), we choose to display also in Table 



XI 



the values for < r{t) > and < s{t) >. With these numbers the reader can readily obtain the values for 
mpcAC with any choice of ca ■ 



C. The J parameter 

In Figures ^ and ^ the vector meson masses and hyperfine splittings are plotted against the corresponding 
pseudoscalar masses for all the datasets. It is difficult to identify an unquenching signal from these plots 
- the data seem to overlay each other. Note that in B, it was reported that there was a tendency for the 
vector mass to increase as the sea quark mass decreases (for fixed pseudoscalar mass). The observations 
for the present matched dataset imply that this may have been due to either an 0{a) effect (since the 
dataset in Q was not fully improved at this level) or a finite volume effect. The conclusion therefore is 
that it is important to run at a fixed a in order to disentangle unquenching effects from lattice artifacts 
or finite volume effects. 

A possible explanation as to why there is no signal of unquenching in our meson spectrum is the following. 
Our matched ensembles are defined to have a common rg value, so any physical quantity that is sensitive 
to this distance scale (and the static quark potential itself) will also, by definition, be matched. Our 
mesons, because they are composed of relatively heavy quarks, are examples of such quantities, and this 
is a possible reason why there is no significant evidence of unquenching effects in the meson spectrum. 

When comparing the experimental data points with the lattice data in Figures H and B we note that the 
lattice data are high. This could be due to an incorrect value of tq being used (rp = 0.49 fm) and that 
the true value of rg is somewhat higher. This possibility is discussed again in the next section. 

A further point regarding hyperfine splitting in Figure g| is that the lattice data for the matched ensembles 
tends to fiatten as the sea quark mass decreases. (The quenched data has a distinctly negative slope, 
whereas the k^''''^ — 0.1350 data is fiat.) Thus the lattice data is tending towards the same behaviour 



as the experimental data which lies on a line with positive slope (independent of the value used for rg). 
This behaviour is apparently spoiled by the unmatched run with k^°^ — 0.1355 (see Figure |6[) which 
has a clear negative slope. However the k^'^^ = 0.1355 data does not satisfy the finite volume bound 



of m (see Sec.VE). One would expect that these finite volume effects would squeeze the vector meson 
state more than the pseudoscalar state (the p is an extended object). Furthermore, as the valence quark 
mass was decreased, the more the vector mass would be raised by finite volume systematics. These 
considerations match with the observed behaviour of the k^'^'^ = 0.1355 data in Figure ||. The JLQCD 
collaboration |29| ] has recently reported on a finite- volume analysis with the same action as used in this 
work. For f3 = 5.2, k^°^ = 0.1350, they found no evidence of finite-volume effects in their 16'^ data 
for either the pseudoscalar or vector meson. It would be interesting to extend this analysis to their 
(3 = 5.2, k'"''' = 0.1355 dataset. 

The J-parameter is defined |4q] as 



, , dMv 

My 



dMlg 



(32) 

K,K' 



In the context of dynamical fermion simulations, this parameter can be calculated in two ways. The 
first is to define a partially quenched J for each value of the sea quark mass. In this case, the derivative 
in (^2|) is with respect to variations in the valence quark mass (with the sea quark mass fixed). The 
second approach is to define J along what we will term the 'unitary' trajectory, i.e. along n^^^ — k™'. In 



Table XII, the results from both methods are given. These values of J are around 25% lower than the 



experimental value Jexpt — 0.48(2)^ 

Finally we note that the physical value of J (i.e. that which most closely follows the procedure used to 
determine the experimental value of Jexpt = 0.48(2)), should be obtained from extrapolating the results 
from the first approach to the physical sea quark masses. We call this the third approach. In order to 
perform this extrapolation, we extrapolate the three matched dynamical J values obtained from the first 
approach linearly in {M^g aiy-j2 ^^ {M^g ary-j2 _ q M^g ^''^ is the pseudoscalar meson mass at the 
unitary point (i.e. where the valence and sea quark masses are all degenerate). The value for J from the 



third approach is presented in Table XII and we note that it is approaching the experimental value for 
J. 

The results from all three approaches are plotted in Figure M, together with the experimental result. 
There is some promising evidence that the lattice estimate of J increases towards the experimental point 
as the sea quark mass decreases (see the J value from approaches 1 and 3). This effect will move the 
lattice estimates of the J parameter towards the experimental value as simulations are performed at more 
physical values of quark mass. 

Recently there has been a proposed ansatz for the functional form of My as a function of Mf^g |47|| . 
However, all our data have Mps/Mv ^ 0.6, and for this region, the ansatz of |43] is linear to a good 
approximation. Therefore we choose to interpolate our data with a simple linear function and await more 
chiral data before using the ansatz of |Q . 

Two groups have recently reported results on the J parameter from dynamical simulations. The CP- 
PACS collaboration results at a w 0.11 fm, found jOynamtcai > jQuenched ^^^^^ ^ ^^^^^^ action |]. 
Furthermore they found that this discrepancy increased as the continuum limit was taken. A similar 
result was found by the MILC collaboration who used an improved staggered action with a ~ 0.13 fm 
iQ. Both these groups' results match those found in this work. 



'Note, however, that the experimental value of J does depend on assumptions regarding the mixing of the 
strange and non-strange quark states. 



D. Lattice spacing 



In Section IV the lattice spacing, a, was determined from the intermediate range properties of the static 
quark potential. In this subsection, we present a complementary determination of a from the meson 
spectrum. 

A common method of determining a from the meson spectrum uses the p mass. However, this requires 
the chiral extrapolation of the vector meson mass down to (almost) the chiral limit. This extrapolation 
is often performed using a linear function. However, as was discussed in the previous subsection, a linear 
chiral extrapolation may not be appropriate for My ^ 0.8 GeV. An alternative method of extracting 
the lattice spacing using the vector meson mass at the simulated data points (i.e. without any chiral 
extrapolation) was given in |48|] . Using this method, we obtain the lattice spacing values as shown in 



Table XIH. Note that these are in general 10-15% larger than the values from Section IV where the 
lattice spacing was determined from vq. A possible explanation for this discrepancy is that the potential 
and mesonic spectrum are contaminated with different 0{a^) errors, or that the value vq = 0.49 fm is 
10-15% too small, and that the true value is tq ss 0.55fm. 

It is interesting to study the lattice spacing determinations in more detail since they are a measure 
of unquenching effects in dynamical simulations. Specifically, it is often assumed that the reason the 
various quenched determinations of a from e.g. the meson spectrum differ from that of rg or the string 
tension is due to dynamical quark effects. An obvious quantity to monitor the merging of the various a 
determinations can be defined: 



aj{/3,rhs 



dij(/9, msca) = 1 - ^ ^a ^ — T' (33) 



where a^ is the lattice spacing determined from the physical quantity i = {Mp, Mk, f-^ . . .}. Obviously, 
if Si J = then the lattice prediction for quantity i using the scale determined from j (or vice versa) is in 
exact agreement with experiment. 

Since our simulations arc improved to 0{a^) we expect that S —>■ 0{a^) as ?nsoa=vai —>■ mi (where mi is 
the average u, d quark mass). Thus a plot of 6 against {Mp'^ aiy-j_2 -^^rould be insightful, where Afp" ^'^ 
is the pseudoscalar mass at the unitary point i.e. for degenerate valence and sea quarks (so Mp^ ^''■^ — oo 
for the quenched data). Here, we work with {M^ ary-j_2 j-g^f^^j- than l/rhsca, for the x— co-ordinate since 
it is equivalent to, but easier to define than l/rhsca- It is important to note that the a;— co-ordinate in 
this plot is the 'control parameter' for the study of unquenching effects, i.e. when we vary this parameter 
from its quenched value towards its experimental value, we hope to see the data plotted in the y— co- 
ordinate move towards its appropriate experimental value. Thus it is easier to interpret unquenching 
effects directly from this plot than from, e.g., plots of My against Mpg for various mgca- 

In Figure H, Si,j is plotted against (Afp" aiy-j_2 ^^^ ^^^^ matched datasets. In this plot we have fixed j = tq 
and the various physical quantities i are ^/a and the hadronic mass pairs {Mk*, Mk) & {Mp, Mt^). The 



method that was used to determine the scale ai from these mass pairs is that of 48 1. It is worth noting 
that the experimental point on this same plot would occur at an a;— co-ordinate of (M^r)"^ ~ 200. 

Figured does not show signs of unquenching for quantities involving the hadronic spectrum i.e. the mass 
pairs {Mk*, Mk) & {Mp, M^^). (Future work will study Si for the matrix element quantities i — f-n and 
/a'.) However, there is evidence of unquenching effects when comparing the scale from tq with that from 
i/a. The quenched value of 5/^ is distinct from the dynamical values, though we note that the method 



used to obtain a was optimised for the extraction of tq rather than a itself (see Sec, IV B 



One may wonder if the the S values may have been distorted by not choosing the simulation parameters 
(/3, rhsea) exactly on the matched trajectory. In order to obtain a rough estimate of the effect of a 
mis- matched value of (3, we use the renormalisation group inspired ansatz for Oj p9| , p0| : 

ar{gl) - ^-'fpAal) X [1 + XJpt(5o')"'] : (34) 



where fprig^) is the usual asymptotic scahng function obtained from integrating the /3-function of QCD 
and Xi is the coefficient of the 0{a") lattice systematic. The functional form for a{gQ) was originally 
applied for the quenched theory, but let us assume that it can also be applied in the unqucnched case. 
Using Eq. (p4) , we see that a mismatch in /3 of A/9 would lead to a relative error in 5 of 

^^^^i||^--3A^. (35) 

This shows that even an error in /3 of as much as A/3 « 0.01 introduces a relative error in S{P) of only 
3%, ruling out any possible mismatching in /? as leading to a significant distortion in S. 



E. Edinburgh Plot 



In Table XIV the ratios Mps/My are displayed for the case k^°^= k™^ The average u and d quark mass 
is fixed by requiring Mps/My = 0.18. As can been seen, the simulations are at much larger dynamical 
quark masses. Figure ^ shows the 'Edinburgh plot' {Mn/Mv v.s. Mps/My) for all the data sets. There 
is no significant variation within the dynamical data as the sea quark mass is changed but the dynamical 
data does tend to lie above the (matched) quenched data. This latter feature may be indicative of finite 
volume effects since these are expected to be larger in full QCD compared to the quenched case [pl|. In 
P] an analysis of dynamical finite volume effects concluded that they were statistically insignificant for 
spatial extents of L ^ 1.6 fm and sea quark masses corresponding to Mps/My ^ 0.67 with around 100 
configurations. This bound is satisfied for the matched ensembles, but not for the k^°^ — 0.1355 case 
where L = 1.60 fm and Mps/My = 0.58. 

Note also that for heavy valence quark masses, the dynamical data lies close to the phenomenological 
curve 153], whereas it tends to drift higher than the curve for small valence quark masses. The (matched) 
quenched data agrees well with the curve. 

Dynamical results for baryons have recently been reported by two groups. CP-PACS (using a clover ac- 
tion) find good agreement with experiment for strange baryons, but their light baryons (in the continuum 
limit) arc around 10% higher than experiment (see sec.VC in j^). They discuss the possibility that this 
is caused by finite volume effects. The MILC collaboration (using an improved staggered action) find 
their dynamical and quenched Edinburgh plots overlay each other p3|. 



F. Chiral Extrapolations 

There are a number of different 'chiral extrapolations' that one can perform in the case of dynamical 
fermions where there is a two-dimensional quark mass parameter space, (rrisca, 'Tivai). We describe three 
such extrapolations of the data. The first uses a partially quenched analysis where each of the nisca 
datasets is extrapolated entirely separately. The second uses only the unitary sub-set with msoa = Tivai. 
The third does a combined fit of all the matched data using a fitting ansatz to model the variation between 
the different rUsea values. 

Note that there have been recent proposals for the functional form of Mjv and My as a function of 
Mfyg which go beyond the usual chiral linear ansatz normally used in extrapolations of lattice data 
|p3,B^. However, as reported in sec.VC, the non-linearity of these functional forms becomes relevant 



only for lattice data lighter than in our simulations and therefore we choose to use naive linear chiral 
extrapolations in the following. 



1. Partially Quenched Chiral Extrapolations 

A partially quenched chiral extrapolation was performed for the hadronic masses M — My , Mn and Ma 
against Mpg, i.e. the following ansatz was used, 

M{P, n'"'^; K™^) = A + BMps{f5, k'"=^; k™^^. (36) 

We have introduced the following nomenclature. In M(/3, k^°^; k™'), the first two arguments refer to the 
sea parameters: the gauge coupling /3 and the sea quark mass k'^'^^. The third argument refers to the 



valence quark mass k™'. The results for these partially quenched extrapolations appear in Table XV. 
Note that there is no convincing sign of unquenching effects in that the A and B values for the matched 
datasets tend to overlay eachother, and there is no clear trend for these values as a function of rusea- 

Although we choose to extrapolate with respect to Mp5(/3, k*"^**; k™')^, we also show, for completeness, 
the values of Kcrit in Table KVl. These were obtained from the usual fit of Afp5'(/3, k^'^''; k™')^ versus 



2. Unitary Chiral Extrapolations 

An extrapolation of the hadronic masses M — My, M^ and Ma against M'jjg was performed for the 
unitary subset of data i.e. the following ansatz was used, 

M(/3, k"'^^; k'^'=^) = A"'"*^'y + B""'*"yMps(/3, k'^''^ te^^f, (37) 

Note that only the matched, dynamical datasets were included in these fits. The results appear in Table 



XVII 



3. Combined Chiral Fits 

It is instructive to perform a combined chiral fit to the entire matched dataset. In order to achieve this 
we consider the following fitting ansatz for the fitting of the hadronic mass M, where 

M(/3, k"°^; k™1) = A'=°™''"'^'* + B'=°™''™^'*Mps(/3, k"""; k™!)^ 

Bo + BiMps{/3, k'^"^ k'"'')-^] MpsiP, k'"''; k™^)^. (38) 

One advantage of such a fitting procedure is that in total, to fit the entire matched dataset, there are 
fewer fitting parameters than are required in the partially quenched analysis. The functional form in 
eq.(^8|) is the simplest functional form which allows for a variation of A and B with the sea quark mass, 
and which is finite for all the datasets studied. (Note that Mps{f3, k^°^; k^'^^) = oo for the quenched data.) 
The other advantage is that the parameters Ai and Bi are a direct measure of unquenching effects. 



The results for the fitting parameters Aq^i and i?o,i are displayed in Table [XVII1| for the hadronic masses 
MvtMm and Ma. The parameters Ai and Bi for all the hadrons are compatible with zero at the 2(j 
level, underlining again the fact that we have not unambiguously uncovered unquenching effects in the 
meson and baryon spectra. 



G. Glueballs and torelons 

Experiment has not so far detected glueball states unambiguously in the hght hadron spectrum. This 
failure is usually believed to be a consequence of mixing between the light glueballs and qq states ( "quarko- 
nia") with the same quantum numbers and similar masses. We lack, however, a clear understanding of 
the mixing matrix elements that lead to the strong interaction cigcnstatcs that would be observed, and 
thus phenomenological attempts to describe the content (gluonic or quarkonium) of the scalar sector 
glueball candidates have led to widely differing results p4| , p5[ . 

Lattice QCD can in principle predict these mixing parameters, and in the quenched approximation 
precise values are known for the continuum gluodynamics (quenched QCD) glueball masses (see ]5^ , |57[ | 
for reviews). Attempts to measure the mixing matrix have been made (see [ p7[ for a review of quenched 
measurements, and [pq] for first determinations in the presence of sea quarks) and arc in progress using the 
current UKQCD field configurations [ p9[ . Simultaneously, the validity of such a simple mixing scenario 
can also be addressed |6Ct| . 

Quenched glueball calculations require large ensembles and, until recently, it had been assumed that a 
similar level of statistical noise would preclude accurate measurements in simulations with dynamical 
fermions. We find, however, and in common with other recent studies [ pl[ that statistical errors are, 
somewhat surprisingly, reduced in dynamical simulation estimates of glueball masses at present parameter 
values, at least compared to similarly sized quenched ensembles. 

Before continuing with a discussion of our calculations we need to be a little more specific about what 
we mean when we talk of "glueballs" in QCD. The point is that the presence of quarks will change the 
vacuum and there is no fundamental reason to think that the mass spectrum of QCD can be approximately 
described as consisting of the glueballs of the pure gauge theory, the usual quarkonia and, where these 
are close in mass, mixtures of the two. There is, however, a collection of phenomena — the OZI rule, 
small sea quark effects etc. — that creates a reasonable prejudice that this might be so. This question 
will be examined more explicitly elsewhere ^^1 ■ Here we shall follow the usual view and assume it to be 
so. In that case we expect that if there are no nearby quarkonia then the states most readily visible using 
purely gluonic operators similar to those used in pure gauge theories will be almost entirely glueball-like. 
This is (probably) the case for the scalar "glueball" state we discuss herein. The fact that the overlap of 
this state onto these purely gluonic operators is similar to that in the pure gauge theory reinforces our 
prejudice. Thus we will refer to this state as the scalar glueball during the remainder of our analysis. 

If we then assume that the glueball spectrum of the dynamical theory is not radically different to that 
of the pure gauge gluodynamics, we expect the lightest states to be the scalar and tensor ground states. 
In terms of the reduced symmetries of the space-time lattice, these correspond to the A\^^ and Tr 



2 



representations of the appropriate cubic group. In the continuum where full rotational symmetry is 
restored, these match onto the J^'~^ — 0++, 2++ states. Given the size of our ensembles, we find it 
difficult to resolve lattice masses much beyond Mq ~ 1.2. In gluodynamics, the heavier tensor state 
has a (continuum extrapolated) mass in units of the Sommer scale around r^Ma — 6. The fo values 
tabulated for our ensembles in Table M thus suggest that the scalar and tensor are the states we are most 
likely to be able to study. 

Using a full arsenal of noise reduction techniques it is now possible to make good estimates of the masses 
of these lightest glueball masses using existing ensembles. In this section we present, as an example, the 
scalar and tensor states extracted from one ensemble, that at (/3, k) — (5.20, 0.13550). Full results for all 
couplings, and giving greater details of methodology, will be reported in pO| . Preliminary results have 
appeared in Q. 

Measurements were made after every tenth HMC trajectory giving an ensemble of 830 configurations, 
which may not be uncorrelated. A jack-knife error analysis was performed using ten bins, each 830 
trajectories in size, which were much larger than the autocorrelation times of the observables. This 



ensured statistically uncorrelated averages for neighbouring bins. 

To reduce statistical errors on mass estimates, operators should have a good overlap onto the ground 
state excitation with the specified quantum numbers. This was achieved in two ways. 

Each operator is based on a traced, closed contour of gauge links, which is gauge invariant. We may 
improve the overlap of these operators onto the ground state excitations by "smearing" and "blocking" 
the links. The former is computationally cheap, but the latter has the advantage of doubling the spatial 
extent of the operator with each iteration. This proves especially useful for measuring wavefunctions that 
are not spherically symmetric, such as the tensor. The details of this procedure will be discussed further 
in@. 

A suite of four glueball operators was constructed in each time-slice of the gauge field configurations 
by summing similarly improved contours in the appropriate symmetry combinations p3l| ]. Overall this 
gave twenty-eight operators per symmetry channel. These were cross correlated and a Ltischer-Wolff 
variational analysis [£2| (for details of the exact procedure see Section 3.2 of |6^) used to extract the 
ground states for each of the lowest momentum combinations of the operators (labelled as P • P = 0, 1 . . . 
where P = p ~ pa). All scalar operators (A^^), for example, were found to have a good overlap 
(typically greater than 0.7) onto the ground state. The robustness of the variational analysis was checked 
by examining the behaviour of individual correlation functions, and of sub-sets of the full operator basis. 
In each case the mass estimates were found to be consistent as expected given the good overlap of all 
operators onto the ground state. 

From correlation functions we may define an effective energy as a function of the Euclidean timelike 
separation t (in lattice units) of the creation and annihilation operators: 

EMt) = -log ^0t(t)O(o)) (2^) 

The effective energies of the non-zero momentum states were converted to effective masses assuming the 
lattice dispersion relation 

EiPf = M^ + J2 sin' (^^) (40) 

The signal from the P • P = 1 channel was found to be particularly useful. The mass of the ground state 
excitation were still small enough for reliable effective energy plateaux to be observed, and statistical 
noise was observed to be only of a similar magnitude to the P • P = channel. For P • P = 2, however, 
the energies of the states were too large to be confidently assessed. Where they could be extracted, they 
showed effective mass plateaux consistent with lower momentum channels. Since they did not improve 
the quality of the fits, however, they were not included. Correlated and uncorrelated plateau fits were 
then carried out using P • P = 0, 1 together. As the former fits differed only within errors, for robustness 
we quote uncorrelated results in this summary. 

In Fig. 11^ we plot the effective masses for the various momentum channels of the scalar glueball. A 
clear plateau is seen in each of the momentum channels. Since these plateaux are compatible, indicating 
a restoration of the continuum Lorentz symmetry, we can combine the lowest momentum channels to 
estimate the pure scalar glueball mass Mq = 0.628(30) in lattice units, or r^Mc = 3.17(15) in units 
of the Sommer scale. We note here that the interpolated quenched glueball mass at this lattice spacing 
is toMg = 3.79 (16) |6y], which is significantly above the scalar mass measured here. There would thus 
appear to be strong evidence for a quenching effect in the scalar glueball channel of QCD. We should 
temper this statement slightly, as there are other possible sources of suppression of the scalar glueball 
mass. Firstly there are finite volume effects which are known to suppress the scalar glueball mass. In 
quenched QCD the principle source of this suppression is the mixing of the glueball with torelon pair 
states, e.g. Q, but we shall demonstrate below that in the present case our lattices are large enough for 
any such effects to be very small. 



More seriously, we do not know the size of this effect in the continuum hmit. In the quenched theory there 
are known to be large scaling violations in the A^^ channel for the Wilson action |6£| with the "scalar 
dip" tending to suppress the mass below the continuum value even at relatively small lattice spacings. 
Without a continuum limit extrapolation of the glueball mass, we cannot here speculate as to the size of 
the corresponding effect in the presence of dynamical fermions, but preliminary work suggests the scalar 
dip may indeed be enhanced in the ensemble considered here m& . 

A similar analysis yields a tensor mass estimate of Mq = 1-28 (9) in lattice units, or tqMg = 6.43 (42) 
in units of the Sommer scale. This is compatible with the interpolated mass in the pure glue theory 
roA/c = 5.91(23). 

Colour flux tubes, analogous to that between a static quark and anti-quark pair, but without source or sink 
can exist on a periodic volume. Rather, the flux tube closes on itself through a spatial boundary (assuming 
it to be in the 'confined' phase), forming what is usually termed a torelon. To a first approximation the 
mass of the lightest such state equals the spatial extent of the lattice multiplied by the energy per unit 
length of the flux tube (the string tension) . In the infinite volume limit such states become very massive 
and decouple from the observed spectrum. 

The vacuum expectation value (VEV) of the Polyakov operator that couples to such a torelon loop is zero 
in the confined phase of gluodynamics, as the loop cannot be broken when no sources in the fundamental 
representation exist. Thus, only a combination of at least two torelons with the appropriate symmetries 
can couple to the particle states in the theory. On lattices small enough that the mass of the lightest 
torelon pair is comparable to the scalar glueball mass we will see significant finite volume effects. 

When light dynamical quarks are present, the torelon becomes unstable to decay. In this case, the 
Polyakov loop operator gains a non-zero expectation value. This is an effect analogous to the string 
breaking seen in the static quark potential measured using Wilson loops and is another explicit signal for 
the presence of light dynamical quarks in these simulations. In addition, it becomes possible for torelon 
states to mix with glueballs. Such states are, of course, lighter than the pairs of torelons that mix in the 
quenched theory, and so we might expect to see finite volume effects on larger lattices in the presence of 
dynamical quarks. 

The Polyakov loop operator is defined as the traced product of links in a line through the periodic spatial 
boundary: 

L 

PM(") = Trnt^p(" + fcA), (41) 

for fjL = 1...3. In order to improve statistics we create a basis of operators using improved spatial links as 
before. 

In Fig. |l^ we plot the vacuum expectation value of the P • P = 0, 1, 2 Polyakov loop operators. From 
momentum conservation we expect the VEVs of the non-zero momentum operators to be zero. This is 
seen to be satisfied within less than two standard deviations in all cases, indicating that the statistical 
errors are under control. It also adds significance to the fact that the P • P = operators have a vacuum 
expectation value that deviates substantially from zero. This is clear evidence of flux tube breaking by 
dynamical quark pair production. 

Fitting effective masses from P-P = 0, 1 after a Liischer-Wolff analysis as before, we estimate the torelon 
mass as Mp = 0.77(5). Including the leading order universal string correction [ |67| , we expect the loop 
mass to vary with the lattice size in D dimensions as 

Mp = ,L-^-^^—^. (42) 

From this we estimate the string tension to be ct = 0.052 (3) or, using the Sommer scale to set physical 
units, y/a — 462 (13) MeV in good agreement with the value quoted in Table 0. 



The mass of the hghtest torelon pair, around twice the torelon mass, is thus clearly too heavy to induce 
finite volume effects. Likewise, finite volume effects from meson exchange through the boundary should 
be small, although we do not consider this process here. The mass of the torelon, on the other hand, is not 
much larger than that of the scalar giucball, and there is a possibility of mixing occurring between the two 
which would lead to a finite volume contamination. We thus perform a variational analysis where we cross 
correlate a basis of eight of the "best" scalar glueball operators with the two "best" torelon operators. 
We find the matrix to be block diagonal within errors, and the two lowest eigenstates match closely the 
original glueball and torelon in mass and operator overlap. Thus this finite volume contamination is 
negligible, something which could have been anticipated from the small size of the Polyakov line VEV. 

In summary, we have presented measurements of the scalar and tensor glueball and torelon masses on 
an ensemble of configurations at (/3, k) = (5.20,0.13550). We find clear signals for the presence of light 
sea quarks, both in a scalar glueball mass that is significantly suppressed below the quenched value at a 
comparable lattice spacing, and in the breaking of the confining flux tube as demonstrated by a non-zero 
expectation value for the spatial Polyakov loop operator. Although non-zero, the smallness of these 
VEVs together with the fact that the torelon and torelon pair masses are significantly larger than the 
scalar glueball mass lead us to believe that the suppression of the scalar glueball mass is not a finite 
volume effect, a conclusion which is reinforced by an explicit mixing analysis. The dependence of these 
effects on the sea quark mass, and whether this effect persists in the continuum limit is not, however, 
resolved here. 



VI. THE TOPOLOGICAL SUSCEPTIBILITY AND f„. 

The ability to access the non-perturbative sectors, and to vary parameters fixed in Nature has made 
lattice Monte Carlo simulation a valuable tool for investigating the role of topological excitations in QCD 
and related theories, and it is these that we now consider. 

In quenched lattice calculations, the continuum topological susceptibility now appears to be relatively free 
of the systematic errors arising from the discretisation, the finite volumes and the various measurement 
algorithms employed. Attempts to measure the microscopic topological structure of the vacuum are also 
well advanced (for a recent review, see p8|). The inclusion of sea quarks in lattice simulations, even at 
the relatively large quark masses currently employed, is numerically extremely expensive, and to avoid 
significant finite volume contamination of the results, the lattice must be relatively coarse, with a spacing 
a ~ 0.1 fm as in this study. Compared to quenched lattice studies at least, this is a significant fraction 
of the mean instanton radius, and has so far precluded a robust, detailed study of the local topological 
features of the vacuum in the presence of sea quarks. The topological susceptibility, on the other hand, 
may be calculated with some confidence and provides one of the first opportunities to test some of the 
more interesting predictions for QCD. Indeed, it is in these measurements that we find some of the most 
striking evidence for the effects of sea quarks (or, alternatively, for a strong quenching effect) in the lattice 
simulations described in this paper. 

We find clear evidence for the expected suppression of the topological susceptibility in the chiral limit, 
despite our relatively large quark masses. From this behaviour we can directly estimate the pion de- 
cay constant without needing to know the lattice operator renormalisation factors that arise in more 
conventional calculations. 

These results were presented at the IOP2000 [Q, the Confinement IV |1^ and, in a much more prelim- 
inary form, the Lattice '99 ||l^] conferences. Since then, we have increased the size of several ensembles 
and included a new parameter set. We also have more accurate results from the quenched theory with 
which to compare. Related results have been presented by the CP-PACS collaboration |69-7l], the Pisa 



group ||7^,[73| , the SESAM-TxL collaboration ||7j and the Boulder group [^ . A detailed analysis of our 
data set, and its relation to these other studies will be given in ||76||. 



The topological charge is 

Q - ^ ^ d^x\e^,arFlMF:,{x). (43) 

The topological susceptibility is the squared expectation value of the topological charge, normalised by 
the volume 

Sea quarks induce an instanton-anti-instanton attraction which in the chiral limit becomes stronger, 
suppressing Q and x [0 

X = sf— + — ) , (45) 

where 

E = - lim lim (O|-0'0|O) (46) 

is the chiral condensate |Q. We assume (0|VJ^|0) = (0|u-u|0) = (0|Jd|0) and neglect contributions of 
heavier quarks. The Gell-Mann-Oakes-Renner relation, 

/2A/2 = 2(m„ + m<i)E + ©(m^) (47) 

implies 

X=^ + 0{M^) (48) 

for Nf degenerate light flavours, in a convention where the experimental value of the pion decay constant 
/^ ~ 132 MeV |. Equation (H) holds in the limit f^M^V > 1, which is satisfied by all our lattices. The 
higher order terms ensure that x ~* X'^'^j the quenched value, as mg,M^ — > cx3. We find, however, that 
our measured values are not very much smaller than x'^", so we must consider two possibilities. 



Firstly, there are phenomenological reasons W9 80| for believing that QCD is 'close' to Nc = 00, and 



in the case of gluodynamics even SU(2) is demonstrably close to SU(oo) [ |64y5qj8l| ]. Fermion effects are 
non-leading in N^, so we expect x ~^ X'^^ for ^W fixed value of niq as the number of colours iVj, — > 00. 
For small mq we expect 

X°°M2 

with x°°, foo the quantities at leading order in N^ [ |7q |. Alternatively, our rUq ~ 'TigtTango and perhaps 
higher order terms arc important. In the absence of a QCD prediction, 

X=^Mjarctan(^X-^) (50) 



** N.B. there is a connnon alternative convention, used in earlier presentations of this data J13| , [l5| , where a factor 
of 2 is absorbed into /^ in mm, and where f^ is a factor of \/2 smaller, around 93 MeV. 



interpolates between (Mq) and the quenched limit .Fj. Measurements of x were made on a number of 
ensembles of Nf = 2 lattice field configurations. We reiterate here that these ensembles have two notable 
features. The improvement is fully non-perturbative, with discretisation errors being quadratic rather 
than linear in the lattice spacing. Second, the couplings are chosen to maintain an approximately constant 
lattice spacing (as defined by the Sommer scale, rg = 0.49 fm [pj) as the quark mass is varied. This is 
important, as the susceptibility in gluodynamics varies considerably with the lattice spacing |56 8^|, in 
competition with the variation with niq . The topological susceptibility is measured from the gauge fields 
after cooling to remove the UV noise. Further details of the procedure may be found in |13|j76|] . 

We plot data for the ensembles presented in this paper in Figs. |l^ and |l^, as well as for preliminary 
results for two further data sets at {(3, n) = (5.20, 0.13565) and (5.25, 0.13520). Also shown, as a band, is 
the interpolated x'^^ ^it an equivalent lattice spacing. Owing to the systematic differences in the methods 
for determining fo (which can amount to a 20% difference in fp), the value chosen is for the quenched 
coupling /3 = 5.93, taken from ||81| , where we have an estimate of fo determined in a consistent manner. 
The variation in the equivalent quenched susceptibility over the range in fo spanned by our data is much 
smaller than the error on the (3 = 5.93 point shown, a useful consequence of the matching programme. 

The behaviour oiff^x with {foAd-^)'^ is qualitatively as expected and, more quantitatively, we attempt fits 
motivated by (Q, (|9|) and (|50|). The leading order chiral behaviour will be 

. 2 - 



CO, (51) 

TT 

with the first correction term generically being 

co + ci(roM,)^ (52) 



Ml 



^O^X _ , / . ,> ^2 



Ml 

Attempting to include data further from the chiral limit, large- iVc theory suggests a functional form 

n)^X ^ C0C3 

Ml c3 + coifoM.r' 

whilst a more general interpolation is provided by 

nj^X 2co / 7rc3 \ 

——TT — arctan . (54) 

Ml TT \2co{foM^)y 

In each case the intercept is related to the decay constant by cq = (fo/7r)^/8. The corresponding fits are 
shown in Figs, g^and |l3|. The extent of the curves indicates which points were included in fit. We include 
progressively less chiral points until the x^/dof of the fit becomes unacceptably bad. We note the wide 
range fitted simply by including an M^ term, and the consistency of our data with large- iVc predictions. 



The stability and similarity of the fits motivates us to use cq from (52) to estimate /^r = 149 ±8 ^14 MeV 



at a lattice spacing a ~ 0.1 fm, with variation between other fits providing the second, systematic error, 
and in good agreement with the experimental value around 132 MeV. 



VII. CONCLUSIONS 

Two particular features distinguish this work from previous published reports on lattice simulations of 
QCD with dynamical fermions. It represents the first presentation of a wide range of results using the 



**Note that, in describing chiral extrapolations, we adopt the common convention of using n to label quantities 
associated with the pseudoscalar channel irrespective of the quark mass 



fully non-perturbativley improved Wilson action. It also demonstrates the value of a new strategy of 
using so-called 'matched ensembles' which allows a more controlled study of unqucnching effects than 
would otherwise be possible at finite lattice spacing. 

We have presented detailed measurements of the static inter-quark potential, light hadron spectrum, 
scalar and tensor glueballs, torelon states and the topological charge and susceptibility. 

From the analysis of these quantities, we have presented significant evidence of effects attributable to 
dynamical effects (two flavours of light quarks) on 



the static inter-quark potential, particularly at short range (section IV C| ) 



the topological susceptibility (section |VlD ; 



We have also seen some evidence of dynamical quark effects in 



• the effective string tension (section VD); 

• the nucleon mass (section |VE|); 



the scalar glueball mass (section VG); 



For the present range of light quark masses (Mt^JMp ^ 0.58) there is no convincing evidence of effects 
on the light meson spectrum. Nor do we see evidence of string breaking, save indirectly in the small, but 
non-zero, VEV of the winding giuonic flux tube (torelon) operator. 

Further analyses of these ensembles and complementary ones being produced by the QCDSF collabora- 



tion [p2 83 1 are underway. 
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/3 


Csw 


#corif. 


^soa 






K^-' 






5.20 


2.0171 


244 


0.13565 


0.13565 










5.20 


2.0171 


832 


0.1355 


0.1355 


0.1350 


0.1345 


0.1340 




5.20 


2.0171 


600 


0.1350 


0.1350 


0.1345 


0.1340 


0.1335 




5.26 


1.9497 


404 


0.1345 


0.1350 


0.1345 


0.1340 


0.1335 




5.29 


1.9192 


404 


0.1340 


0.1350 


0.1345 


0.1340 


0.1335 




5.93 


1.82 


623 





0.1339 


0.1337 


0.1334 


0.1332 


0.1327 



TABLE I. Summary of simulation parameters and statistics used in the computation of the static potential 
and light hadron spectrum. 



#traj. r'"^ 

16^ • 32 5.20 2.0171 0.13565 2400 13(5) 

14(1) 
16(3) 
18(3) 
25(7) 



5.20 


2.0171 


0.13565 


2400 


5.20 


2.0171 


0.13550 


8000 


5.20 


2.0171 


0.13500 


6000 


5.26 


1.9497 


0.13450 


6000 


5.29 


1.9192 


0.13400 


5000 



16^ • 24 


5.20 


1.76 


0.1390 


3800 


37(3) 




5.20 


1.76 


0.1395 


3200 


27(18) 




5.20 


1.76 


0.1398 


3000 


32(8) 



TABLE II. Comparison of integrated autocorrelation times r'"' for the average plaquette measured in the 
present simulations with those in previous simulations at /3 = 5.20, Csw = 1.76. 



(13, n) L/ro LM^ 

(5.20, 0.13565) 3.07 (3) 4.18 (5) 

(5.20, 0.13550) 3.17 (3) 4.70 (6) 

(5.20, 0.13500) 3.37 (3) 6.48 (8) 

(5.26, 0.13450) 3.40 (4) 8.14 (3) 

(5.29, 0.13400) 3.32 (3) 9.23 (4) 

TABLE III. Measures of finite volume effects in simulations. 



{13, ^=-) 


r 


(5.20, .1350) 


(5.26, .1345 


) (5.29, .1340) 


(5.93,0) 


(5.2, .1355) 


(5.2, .13565) 


(1,0,0) 


0.4823(02) 


0.4739(04) 


0.4707(03) 


0.4259(01) 


0.4762(02) 


0.4749(02) 


(2,0,0) 


0.6970(08) 


0.6839(11) 


0.6782(10) 


0.6268(03) 


0.6832(08) 


0.6794(06) 


(3,0,0) 


0.8253(17) 


0.8100(15) 


0.8027(17) 


0.7439(05) 


0.7999(14) 


0.7954(12) 


(4,0,0) 


0.9193(22) 


0.9001(27) 


0.8920(28) 


0.8307(06) 


0.8839(18) 


0.8745(14) 


(5,0,0) 


0.9945(30) 


0.9777(36) 


0.9654(36) 


0.9070(07) 


0.9504(28) 


0.939(02) 


(6,0,0) 


1.0628(43) 


1.042(06) 


1.0342(43) 


0.9780(09) 


1.0168(29) 


1.002(02) 


(7,0,0) 


1.130(06) 


1.105(06) 


1.098(07) 


1.0484(13) 


1.0828(39) 


1.061(04) 


(8,0,0) 


1.183(08) 


1.175(09) 


1.170(11) 


1.1117(16) 


1.135(05) 


1.114(04) 


(9,0,0) 


1.262(11) 


1.244(11) 


1.244(11) 


1.1802(26) 


1.186(07) 


1.165(05) 


(10,0,0) 


1.321(17) 


1.285(21) 


1.310(15) 


1.243(4) 


1.246(08) 


1.221(07) 


(11,0,0) 


1.398(21) 


1.414(23) 


1.367(16) 


1.301(5) 


1.298(10) 


1.277(11) 


(12,0,0) 


1.467(24) 


- 


- 


1.365(8) 


1.330(17) 


1.287(25) 


(1,1,0) 


0.6276(05) 


0.6156(06) 


0.6103(07) 


0.5514(2) 


0.6173(05) 


0.6140(04) 


(2,1,0) 


0.7495(09) 


0.7315(13) 


0.7288(11) 


0.6671(4) 


0.7310(09) 


0.7262(07) 


(2,2,0) 


0.8163(14) 


0.8001(17) 


0.7940(15) 


0.7319(5) 


0.7944(10) 


0.7884(10) 


(3,1,0) 


0.8483(15) 


0.8296(16) 


0.8226(16) 


0.7616(6) 


0.8215(15) 


0.8138(11) 


(3,2,0) 


0.8873(18) 


0.8687(27) 


0.8636(23) 


0.8009(7) 


0.8599(15) 


0.8497(14) 


(3,3,0) 


0.9387(24) 


0.9235(26) 


0.9122(22) 


0.8517(9) 


0.9051(17) 


0.8939(18) 


# conf. 


150 


101 


101 


623 


208 


244 


traj. spac. 


40 


40 


40 


- 


40 


10x2 


TABLE IV. 


The static potential V(r) in lattice units. For the 


preliminary data at k*"^^ = 0.13565 the configu- 


rations were measured every 10 


trajectories ; 


md analysed in bins 


of 2. 






(/3,^=-) 


ro/a 


a 


[/H 


e 




^[MeV] 


(5.2, .13565) 


5.21(05)(+0- 


8) 0.0941(8)(+13-0) 


0.315(7)(+18-11) 


465(l)(+19-3) 


(5.2, .13550) 


5.041(40)(+0 - 


-10) 


.0972(8)(+7-0) 


0.307(6)(+17-l) 


467(l)(+17-3) 


(5.20, .1350) 


4.754(40)(+2 - 


-90) 0.1031(09)(+20 - 1) 


0.326(07)(+32-12) 


463(2)(+2-6) 


(5.26, .1345) 


4.708(52)(+45 


-50) 0.1041(12)(+11 - 10) 


0.298(09)(+100-8) 


468(2)(+2-18) 


(5.29, .1340) 


4.813(45)(+35 


-84) 


.1018(10)(+20-7) 


0.310(10)(+0-61) 


466(2)(+10-0) 


(5.93,0) 


4.714(13)(+0- 


-18) 


.1040(03)(+4-0) 


0.276(03)(+17-2) 


471(l)(+21-3) 



TABLE V. Sommer scale ro and other parameters deduced from the lattice potential. 



roMp 



aMp 



5.2000 
5.2000 
5.2000 
5.2000 
5.2000 
5.2000 
5.2000 
5.2000 
5.2000 
5.2000 



0.1355 
0.1355 
0.1355 
0.1355 
0.1355 
0.1355 
0.1355 
0.1355 
0.1355 
0.1355 



0.1340 
0.1345 
0.1345 
0.1350 
0.1350 
0.1350 
0.1355 
0.1355 
0.1355 
0.1355 



0.1340 
0.1340 
0.1345 
0.1340 
0.1345 
0.1350 
0.1340 
0.1345 
0.1350 
0.1355 



2.39 + 2 

2.25 + 2 

2.12 + I 

2.12 + 2 

1.97 + 2 
1.82 + I 

1.98 + I 
1.82 + I 
1.66 + I 
1.48 + I 



0.473 
0.447 
0.420 
0.420 
0.391 
0.362 
0.392 
0.362 
0.329 
0.294 



+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 3 

- 2 
+ 3 

- 3 
+ 3 

- 2 
+ 3 

- 3 
+ 3 

- 3 
+ 4 

- 3 



5.2000 


0.1350 


0.1335 


0.1335 


2.68 


5.2000 


0.1350 


0.1340 


0.1335 


2.56 


5.2000 


0.1350 


0.1340 


0.1340 


2.45 


5.2000 


0.1350 


0.1345 


0.1335 


2.44 


5.2000 


0.1350 


0.1345 


0.1340 


2.32 


5.2000 


0.1350 


0.1345 


0.1345 


2.20 


5.2000 


0.1350 


0.1350 


0.1335 


2.32 


5.2000 


0.1350 


0.1350 


0.1340 


2.20 


5.2000 


0.1350 


0.1350 


0.1345 


2.06 


5.2000 


0.1350 


0.1350 


0.1350 


1.93 



+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 



0.563 
0.539 
0.514 
0.514 
0.489 
0.462 
0.488 
0.462 
0.434 
0.405 



+ 3 

- 3 

+ 3 

- 4 
+ 3 

- 4 
+ 3 

- 4 
+ 3 

- 4 
+ 4 

- 5 
+ 3 

- 4 
+ 4 

- 5 
+ 4 

- 5 
+ 4 

- 5 



5.2600 
5.2600 
5.2600 
5.2600 
5.2600 
5.2600 
5.2600 
5.2600 
5.2600 
5.2600 



0.1345 
0.1345 
0.1345 
0.1345 
0.1345 
0.1345 
0.1345 
0.1345 
0.1345 
0.1345 



0.1335 
0.1340 
0.1340 
0.1345 
0.1345 
0.1345 
0.1350 
0.1350 
0.1350 
0.1350 



0.1335 
0.1335 
0.1340 
0.1335 
0.1340 
0.1345 
0.1335 
0.1340 
0.1345 
0.1350 



2.85 + I 

2.74 + I 

2.63 + I 

2.63 + I 

2.52 + I 

2.41 + I 

2.52 + I 

2.41 + I 

2.29 + I 

2.16 + I 



0.603 
0.580 
0.557 
0.557 
0.533 
0.509 
0.533 
0.509 
0.484 
0.458 



+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 
+ 2 

- 2 



5.2900 



0.1340 



0.1335 



0.1335 



2.99 



0.621 



5.2900 0.1340 0.1340 0.1335 



2.88 + 1 0.599 + 2 



5.2900 0.1340 0.1340 0.1340 2.78 + I 0.577 



+ 2 



- 4 "■'-" ' - 2 

5.2900 0.1340 0.1345 0.1335 2.78 + I 0.577 + I 

5.2900 0.1340 0.1345 0.1340 2.67 + I 0.554 + I 

5.2900 0.1340 0.1345 0.1345 2.55 + I 0.530 + 3 

2.67 + I 0.554 + 2 



5.2900 0.1340 0.1350 0.1335 

5.2900 0.1340 0.1350 0.1340 

5.2900 0.1340 0.1350 0.1345 

5.2900 0.1340 0.1350 0.1350 2.31 + ^ 0.480 + ^ 



2.55 + I 0.530 + 3 

2.43 + I 0.506 + 3 

3 0-480 + I 



5.9300 


0.0000 


0.1327 


0.1327 


2.334 


+ 


6 
10 


0.495 


+ 1 
- 1 


5.9300 


0.0000 


0.1332 


0.1327 


2.211 


+ 


6 
9 


0.469 


+ 1 
- 1 


5.9300 


0.0000 


0.1332 


0.1332 


2.081 


+ 


6 
9 


0.442 


+ 1 
- 1 


5.9300 


0.0000 


0.1334 


0.1327 


2.159 


+ 


6 
9 


0.458 


+ 1 
- 1 


5.9300 


0.0000 


0.1334 


0.1332 


2.028 


+ 


6 
9 


0.430 


+ 1 
- 1 


5.9300 


0.0000 


0.1334 


0.1334 


1.973 


+ 


6 
9 


0.419 


+ 1 
- 1 


5.9300 


0.0000 


0.1337 


0.1337 


1.800 


+ 


6 
9 


0.382 


+ 1 
- 1 


5.9300 


0.0000 


0.1339 


0.1337 


1.739 


+ 


6 
9 


0.369 


+ 1 
- 1 


5.9300 


0.0000 


0.1339 


0.1339 


1.676 


+ 


6 
9 


0.356 


+ 1 
- 1 



TABLE VI. Pseudoscalar meson masses for all data sets. 



/3 


^sea 




^val 


TqMv 




aM 


V 


5.2000 


0.1355 


0.1340 


0.1340 


3.01 


+ 


5 
2 


0.596 


+ 6 
- 5 


5.2000 


0.1355 


0.1345 


0.1340 


2.92 


+ 


5 
3 


0.578 


+ 6 
- 6 


5.2000 


0.1355 


0.1345 


0.1345 


2.82 


+ 


5 
3 


0.560 


+ 7 
- 6 


5.2000 


0.1355 


0.1350 


0.1340 


2.84 


+ 


5 
3 


0.563 


+ 7 
- 6 


5.2000 


0.1355 


0.1350 


0.1345 


2.75 


+ 


6 
3 


0.546 


+ 8 

— 7 


5.2000 


0.1355 


0.1350 


0.1350 


2.68 


+ 


7 
4 


0.531 


+ 10 

- 8 


5.2000 


0.1355 


0.1355 


0.1340 


2.79 


+ 


6 

4 


0.553 


+ 10 

- 8 


5.2000 


0.1355 


0.1355 


0.1345 


2.71 


+ 


7 
4 


0.537 


+ 11 
- 9 


5.2000 


0.1355 


0.1355 


0.1350 


2.63 


+ 


8 
5 


0.522 


+ 13 
- 9 


5.2000 


0.1355 


0.1355 


0.1355 


2.56 


+ 10 

- 4 


0.508 


+ 18 
-10 


5.2000 


0.1350 


0.1335 


0.1335 


3.31 


+ 


3 

4 


0.695 


+ 4 
- 4 


5.2000 


0.1350 


0.1340 


0.1335 


3.22 


+ 


3 

4 


0.677 


+ 4 
- 5 


5.2000 


0.1350 


0.1340 


0.1340 


3.13 


+ 


3 

4 


0.658 


+ 5 
- 5 


5.2000 


0.1350 


0.1345 


0.1335 


3.13 


+ 


3 

4 


0.658 


+ 5 
- 5 


5.2000 


0.1350 


0.1345 


0.1340 


3.04 


+ 


3 

4 


0.638 


+ 5 
- 6 


5.2000 


0.1350 


0.1345 


0.1345 


2.94 


+ 


3 

4 


0.619 


+ 6 

— 7 


5.2000 


0.1350 


0.1350 


0.1335 


3.03 


+ 


3 

4 


0.638 


+ 5 
- 5 


5.2000 


0.1350 


0.1350 


0.1340 


2.94 


+ 


3 

4 


0.618 


+ 6 
- 6 


5.2000 


0.1350 


0.1350 


0.1345 


2.85 


+ 


3 
5 


0.599 


+ 6 

- 7 


5.2000 


0.1350 


0.1350 


0.1350 


2.75 


+ 


4 
5 


0.579 


+ 7 
- 9 


5.2600 


0.1345 


0.1335 


0.1335 


3.41 


+ 


3 
5 


0.721 


+ 4 
- 4 


5.2600 


0.1345 


0.1340 


0.1335 


3.32 


+ 


3 
5 


0.703 


+ 4 
- 4 


5.2600 


0.1345 


0.1340 


0.1340 


3.24 


+ 


3 
5 


0.685 


+ 4 
- 4 


5.2600 


0.1345 


0.1345 


0.1335 


3.24 


+ 


3 
5 


0.685 


+ 4 
- 4 


5.2600 


0.1345 


0.1345 


0.1340 


3.16 


+ 


3 
5 


0.668 


+ 4 
- 4 


5.2600 


0.1345 


0.1345 


0.1345 


3.07 


+ 


3 
5 


0.650 


+ 4 
- 4 


5.2600 


0.1345 


0.1350 


0.1335 


3.16 


+ 


3 
5 


0.668 


+ 4 
- 5 


5.2600 


0.1345 


0.1350 


0.1340 


3.08 


+ 


3 
5 


0.651 


+ 4 
- 5 


5.2600 


0.1345 


0.1350 


0.1345 


2.99 


+ 


3 
5 


0.633 


+ 5 
- 5 


5.2600 


0.1345 


0.1350 


0.1350 


2.91 


+ 


3 
5 


0.614 


+ 5 
- 5 



5.2900 



0.1340 



0.1335 



0.1335 



3.49 



0.725 



5.2900 0.1340 0.1340 0.1335 3.41 + l 0.708 + ^ 



5.2900 0.1340 0.1340 0.1340 3.32 + l 0.691 + « 

— 6 — 6 

5.2900 0.1340 0.1345 0.1335 3.32 + l 0.691 + I 

— 6 — 6 

5.2900 0.1340 0.1345 0.1340 3.24 ^ t 0.674 + I 

— 6 — 6 

5.2900 0.1340 0.1345 0.1345 3.16 + g 0.656 + ^ 

5.2900 0.1340 0.1350 0.1335 3.24 + g 0.674 + ^ 

5.2900 0.1340 0.1350 0.1340 3.16 + t 0.656 + I 

— 6 — 7 

5.2900 0.1340 0.1350 0.1345 3.08 t g 0.639 + I 

5.2900 0.1340 0.1350 0.1350 3.00 + f 0.623 + I 



5.9300 


0.0000 


0.1327 


0.1327 


3.05 


+ 


1 

2 


0.646 


+ 2 
- 3 


5.9300 


0.0000 


0.1332 


0.1327 


2.97 


+ 


1 

2 


0.629 


+ 3 
- 3 


5.9300 


0.0000 


0.1332 


0.1332 


2.88 


+ 


2 
2 


0.612 


+ 3 
- 3 


5.9300 


0.0000 


0.1334 


0.1327 


2.93 


+ 


2 
2 


0.622 


+ 3 
- 3 


5.9300 


0.0000 


0.1334 


0.1332 


2.85 


+ 


2 
2 


0.605 


+ 3 
- 3 


5.9300 


0.0000 


0.1334 


0.1334 


2.82 


+ 


2 
2 


0.598 


+ 3 
- 4 


5.9300 


0.0000 


0.1337 


0.1337 


2.72 


+ 


2 
2 


0.577 


+ 4 
- 4 


5.9300 


0.0000 


0.1339 


0.1337 


2.69 


+ 


2 
2 


0.570 


+ 5 
- 4 


5.9300 


0.0000 


0.1339 


0.1339 


2.66 


+ 


3 
3 


0.563 


+ 5 
- 5 



TABLE VII. Vector meson masses for all data sets. 



/3 




K-"''- 


roMiv 




aM 


'iV 


5.2000 


0.1355 


0.1340 


4.75 


+ 


9 
6 


0.942 


+ 12 
-13 


5.2000 


0.1355 


0.1345 


4.42 


+ 


9 
6 


0.876 


+ 15 
-15 


5.2000 


0.1355 


0.1350 


4.09 


+ 10 

- 7 


0.81 


+ 2 
- 2 


5.2000 


0.1355 


0.1355 


3.86 


+ 


7 
5 


0.766 


+ 10 
-11 


5.2000 


0.1350 


0.1335 


5.16 


+ 


5 
6 


1.086 


+ 8 
- 8 


5.2000 


0.1350 


0.1340 


4.87 


+ 


5 
6 


1.024 


+ 8 
- 9 


5.2000 


0.1350 


0.1345 


4.54 


+ 


5 

7 


0.954 


+ 8 
-11 


5.2000 


0.1350 


0.1350 


4.20 


+ 


5 
7 


0.883 


+ 10 
-12 


5.2600 


0.1345 


0.1335 


5.32 


+ 


5 
9 


1.125 


+ 8 
- 8 


5.2600 


0.1345 


0.1340 


5.05 


+ 


5 
9 


1.068 


+ 9 

- 8 


5.2600 


0.1345 


0.1345 


4.78 


+ 


5 
9 


1.011 


+ 10 
- 9 


5.2600 


0.1345 


0.1350 


4.50 


+ 


6 
9 


0.951 


+ 10 
-10 


5.2900 


0.1340 


0.1335 


5.50 


+ 


5 
9 


1.143 


+ 8 
- 8 


5.2900 


0.1340 


0.1340 


5.23 


+ 


5 
8 


1.086 


+ 9 
- 9 


5.2900 


0.1340 


0.1345 


4.94 


+ 


6 
9 


1.027 


+ 10 
-10 


5.2900 


0.1340 


0.1350 


4.66 


+ 


7 
9 


0.968 


+ 13 
-12 


5.9300 


0.0000 


0.1327 


4.56 


+ 


2 
3 


0.968 


+ 5 
- 6 


5.9300 


0.0000 


0.1332 


4.25 


+ 


3 

4 


0.902 


+ 5 
- 8 


5.9300 


0.0000 


0.1334 


4.13 


+ 


3 

4 


0.876 


+ 6 

- 8 


5.9300 


0.0000 


0.1337 


3.94 


+ 


3 

5 


0.836 


+ 7 
- 9 


5.9300 


0.0000 


0.1339 


3.86 


+ 


4 
4 


0.818 


+ 7 
- 8 



TABLE VIII. Nucleon masses for all data sets. 



/3 




K-"''- 


roMA 


uMa 


5.2000 


0.1355 


0.1340 


5.12 


+ 10 

- 5 


1.015 


+ 15 
-12 


5.2000 


0.1355 


0.1345 


4.87 


+ 10 
- 6 


0.967 


+ 17 
-15 


5.2000 


0.1355 


0.1350 


4.64 


+ 11 

- 7 


0.92 


+ 2 
- 2 


5.2000 


0.1355 


0.1355 


4.30 


+ 15 
-11 


0.85 


+ 3 
- 2 


5.2000 


0.1350 


0.1335 


5.57 


+ 6 
- 8 


1.172 


+ 11 
-11 


5.2000 


0.1350 


0.1340 


5.31 


+ 6 

- 8 


1.116 


+ 11 
-12 


5.2000 


0.1350 


0.1345 


5.02 


+ 7 
- 8 


1.055 


+ 13 

-15 


5.2000 


0.1350 


0.1350 


4.75 


+ 8 
-10 


1.00 


+ 2 
- 2 


5.2600 


0.1345 


0.1335 


5.61 


+ 5 
- 9 


1.186 


+ 11 
-10 


5.2600 


0.1345 


0.1340 


5.36 


+ 5 
- 9 


1.134 


+ 11 
-11 


5.2600 


0.1345 


0.1345 


5.11 


+ 6 
- 9 


1.080 


+ 12 
-11 


5.2600 


0.1345 


0.1350 


4.83 


+ 7 
-10 


1.022 


+ 14 
-13 


5.2900 


0.1340 


0.1335 


5.80 


+ 6 
-10 


1.205 


+ 11 
-10 


5.2900 


0.1340 


0.1340 


5.56 


+ 6 
-10 


1.155 


+ 12 
-11 


5.2900 


0.1340 


0.1345 


5.33 


+ 6 
-10 


1.107 


+ 11 
-13 


5.2900 


0.1340 


0.1350 


5.09 


+ 7 
- 9 


1.057 


+ 13 
-12 


5.9300 


0.0000 


0.1327 


5.09 


+ 3 

- 4 


1.079 


+ 7 
- 8 


5.9300 


0.0000 


0.1332 


4.84 


+ 4 
- 5 


1.026 


+ 8 
- 9 


5.9300 


0.0000 


0.1334 


4.74 


+ 4 
- 5 


1.005 


+ 9 
- 9 


5.9300 


0.0000 


0.1337 


4.58 


+ 5 
- 5 


0.972 


+ 11 
-11 


5.9300 


0.0000 


0.1339 


4.47 


+ 6 
- 6 


0.949 


+ 12 
-11 



TABLE IX. Delta masses for all data sets. 



rompcAC ampcAC 



5.20 


0.1355 


0.1340 


0.1340 


0.329 


5.20 


0.1355 


0.1345 


0.1340 


0.292 


5.20 


0.1355 


0.1345 


0.1345 


0.256 


5.20 


0.1355 


0.1350 


0.1340 


0.256 


5.20 


0.1355 


0.1350 


0.1345 


0.221 


5.20 


0.1355 


0.1350 


0.1350 


0.185 


5.20 


0.1355 


0.1355 


0.1340 


0.221 


5.20 


0.1355 


0.1355 


0.1345 


0.186 


5.20 


0.1355 


0.1355 


0.1350 


0.151 


5.20 


0.1355 


0.1355 


0.1355 


0.116 



+ 2 0.0652 + I 

+ 2 0.0580 + 3 

+ 2 0.0508 + 2 

2 0.0508 + 3 

2 0.0438 + 3 

0.0368 + 3 

+ ^ 0.0438 + 3 



- 2 
+ 3 

- 2 



_ 2 u.u^oo _ 3 

+ 3 

- 2 
+ 2 

- 2 
+ 2 

- 2 



0.0368 + 3 
0.0299 + 3 
0.0231 + 3 



5.20 


0.1350 


0.1335 


0.1335 


0.424 


5.20 


0.1350 


0.1340 


0.1335 


0.389 


5.20 


0.1350 


0.1340 


0.1340 


0.355 


5.20 


0.1350 


0.1345 


0.1335 


0.355 


5.20 


0.1350 


0.1345 


0.1340 


0.320 


5.20 


0.1350 


0.1345 


0.1345 


0.287 


5.20 


0.1350 


0.1350 


0.1335 


0.320 


5.20 


0.1350 


0.1350 


0.1340 


0.286 


5.20 


0.1350 


0.1350 


0.1345 


0.253 


5.20 


0.1350 


0.1350 


0.1350 


0.220 



+ 3 
4 



0.0893 + 2 

+ I 0.0819 + 2 

+ I 0.0746 + 2 

+ ^ 0.0746 + 2 



_ 4 u.ui^u _ T 

+ 3 0.0674 + 2 



- 3 

+ 2 

- 2 



0.0602 + 3 

3 0.0674 + 2 

3 0.0602 + 3 

0.0532 + 3 

0.0462 + 3 



- 3 u.i^uuz. _ 3 

+ ^ 0.0532 + 3 



5.26 


0.1345 


0.1335 


0.1335 


0.491 


5.26 


0.1345 


0.1340 


0.1335 


0.455 


5.26 


0.1345 


0.1340 


0.1340 


0.420 


5.26 


0.1345 


0.1345 


0.1335 


0.420 


5.26 


0.1345 


0.1345 


0.1340 


0.385 


5.26 


0.1345 


0.1345 


0.1345 


0.351 


5.26 


0.1345 


0.1350 


0.1335 


0.385 


5.26 


0.1345 


0.1350 


0.1340 


0.351 


5.26 


0.1345 


0.1350 


0.1345 


0.317 


5.26 


0.1345 


0.1350 


0.1350 


0.283 



I 0.1038 + 3 

0.0963 + 3 

0.0888 + 3 



_ Q u.iuoo _ 3 

+ 3 

- 7 ^■^■^^^ _ 3 

+ 3 

- 6 

+ I 0.0888 + I 

- 6 — 3 

+ g 0.0815 + 3 



+ 3 

- 6 

+ I 0.0814 + I 

— 6 — 3 

+ 3 

- 6 

+ I 0.0670 + I 

— 5 — 3 



0.0742 + 3 



0.0742 + 3 



0.0599 + 3 



5.29 0.1340 0.1335 0.1335 0.530 + t 0.1101 + ^ 



7 ^.-.-.v^-. _ 3 



5.29 0.1340 0.1340 0.1335 0.494 + ^ 0.1026 + 3 

5.29 0.1340 0.1340 0.1340 0.458 + t 0.0952 + i? 

— 6 — 3 

5.29 0.1340 0.1345 0.1335 0.458 + t 0.0951 + I 

— o —3 

5.29 0.1340 0.1345 0.1340 0.423 + t 0.0878 + I 

— 6 — 3 

5.29 0.1340 0.1345 0.1345 0.387 + I 0.0805 + I 

— 5 — 3 

5.29 0.1340 0.1350 0.1335 0.422 + t, 0.0877 + I 

— 6 — 3 

5.29 0.1340 0.1350 0.1340 0.387 ■*" I 0.0805 + I 

— 5 — 3 

5.29 0.1340 0.1350 0.1345 0.353 "*" I 0.0733 + i? 

— 5 — 3 

5.29 0.1340 0.1350 0.1350 0.318 + I 0.0661 + I 

— 5 — 3 



5.93 


0.0000 


0.1327 


0.1327 


0.3530 


+ 10 
-12 


0.07488 


+ 11 
-11 


5.93 


0.0000 


0.1332 


0.1327 


0.3162 


+ 9 
-11 


0.06709 


+ 11 
-11 


5.93 


0.0000 


0.1332 


0.1332 


0.2799 


+ 8 
-10 


0.05938 


+ 11 
-11 


5.93 


0.0000 


0.1334 


0.1327 


0.3016 


+ 9 
-11 


0.06398 


+ 11 
-12 


5.93 


0.0000 


0.1334 


0.1332 


0.2653 


+ 8 
-10 


0.05629 


+ 11 
— 12 


5.93 


0.0000 


0.1334 


0.1334 


0.2508 


+ 8 
-10 


0.05322 


+ 12 
-11 


5.93 


0.0000 


0.1337 


0.1337 


0.2077 


+ 7 
- 8 


0.04406 


+ 12 
-11 


5.93 


0.0000 


0.1339 


0.1337 


0.1931 


+ 7 
- 8 


0.04097 


+ 12 
-11 


5.93 


0.0000 


0.1339 


0.1339 


0.1786 


+ 7 
- 8 


0.03788 


+ 13 
-12 



TABLE X. The quark mass mpcAC as defined in eq.(^0|) for all the datasets. 



r{t) 



sit) 



5.20 


0.1355 


0.1340 


0.1340 


0.0662 


5.20 


0.1355 


0.1345 


0.1340 


0.0589 


5.20 


0.1355 


0.1345 


0.1345 


0.0516 


5.20 


0.1355 


0.1350 


0.1340 


0.0516 


5.20 


0.1355 


0.1350 


0.1345 


0.0445 


5.20 


0.1355 


0.1350 


0.1350 


0.0374 


5.20 


0.1355 


0.1355 


0.1340 


0.0445 


5.20 


0.1355 


0.1355 


0.1345 


0.0374 


5.20 


0.1355 


0.1355 


0.1350 


0.0304 


5.20 


0.1355 


0.1355 


0.1355 


0.0235 



+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 2 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 2 

- 3 
+ 3 

- 3 
+ 3 

- 3 



0.1179 


+ 9 

- 8 


0.1050 


+ 9 

- 8 


0.0923 


+ 9 

- 8 


0.0924 


+ 9 
- 8 


0.0799 


+ 10 
- 8 


0.0678 


+ 9 
- 8 


0.0801 


+ 9 
- 8 


0.0678 


+ 9 

- 8 


0.0558 


+ 9 

- 8 


0.0441 


+ 9 



5.20 


0.1350 


0.1335 


0.1335 


0.0907 


5.20 


0.1350 


0.1340 


0.1335 


0.0832 


5.20 


0.1350 


0.1340 


0.1340 


0.0758 


5.20 


0.1350 


0.1345 


0.1335 


0.0758 


5.20 


0.1350 


0.1345 


0.1340 


0.0685 


5.20 


0.1350 


0.1345 


0.1345 


0.0612 


5.20 


0.1350 


0.1350 


0.1335 


0.0685 


5.20 


0.1350 


0.1350 


0.1340 


0.0612 


5.20 


0.1350 


0.1350 


0.1345 


0.0541 


5.20 


0.1350 


0.1350 


0.1350 


0.0469 



0.1682 
0.1538 
0.1397 
0.1397 
0.1258 
0.1123 
0.1260 
0.1123 
0.0990 
0.0859 



+ 11 

- 9 
+ 10 

- 9 
+ 10 
-10 
+ 10 
-10 
+ 10 

- 9 
+ 10 

- 9 
+ 10 
-10 
+ 10 
-10 
+ 10 
-10 
+ 10 

- 9 



5.26 


0.1345 


0.1335 


0.1335 


0.1055 


5.26 


0.1345 


0.1340 


0.1335 


0.0978 


5.26 


0.1345 


0.1340 


0.1340 


0.0903 


5.26 


0.1345 


0.1345 


0.1335 


0.0902 


5.26 


0.1345 


0.1345 


0.1340 


0.0828 


5.26 


0.1345 


0.1345 


0.1345 


0.0754 


5.26 


0.1345 


0.1350 


0.1335 


0.0827 


5.26 


0.1345 


0.1350 


0.1340 


0.0753 


5.26 


0.1345 


0.1350 


0.1345 


0.0680 


5.26 


0.1345 


0.1350 


0.1350 


0.0608 



+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 
+ 3 

- 3 



0.1924 
0.1779 
0.1636 
0.1637 
0.1496 
0.1359 
0.1498 
0.1360 
0.1225 
0.1093 



+ 11 

- 9 

+ 10 

- 9 
+ 10 

- 8 
+ 10 

- 8 
+ 9 

- 8 
+ 9 

- 8 
+ 10 

- 8 
+ 9 

- 8 
+ 9 

- 8 
+ 9 



5.29 



0.1340 



0.1335 



0.1335 



0.1119 



0.2005 



5.29 


0.1340 


0.1340 


0.1335 


0.1042 


+ 3 
- 3 


0.1862 


+ 12 
-12 


5.29 


0.1340 


0.1340 


0.1340 


0.0967 


+ 3 
- 3 


0.1721 


+ 12 
-11 


5.29 


0.1340 


0.1345 


0.1335 


0.0966 


+ 3 
- 3 


0.1721 


+ 12 
-12 


5.29 


0.1340 


0.1345 


0.1340 


0.0892 


+ 3 
- 3 


0.1583 


+ 12 
-12 


5.29 


0.1340 


0.1345 


0.1345 


0.0818 


+ 3 
- 3 


0.1447 


+ 12 
-11 


5.29 


0.1340 


0.1350 


0.1335 


0.0891 


+ 3 
- 3 


0.1583 


+ 12 
-12 


5.29 


0.1340 


0.1350 


0.1340 


0.0817 


+ 3 
- 3 


0.1447 


+ 12 
-12 


5.29 


0.1340 


0.1350 


0.1345 


0.0744 


+ 3 
- 3 


0.1314 


+ 12 
-11 


5.29 


0.1340 


0.1350 


0.1350 


0.0671 


+ 3 
- 3 


0.1182 


+ 12 
-11 


5.93 


0.0000 


0.1327 


0.1327 


0.07584 


+ 11 
-12 


0.1260 


+ 4 
- 5 


5.93 


0.0000 


0.1332 


0.1327 


0.06795 


+ 12 
-12 


0.1127 


+ 4 
- 5 


5.93 


0.0000 


0.1332 


0.1332 


0.06014 


+ 11 
-12 


0.0997 


+ 4 
- 5 


5.93 


0.0000 


0.1334 


0.1327 


0.06480 


+ 11 
-12 


0.1075 


+ 4 
- 5 


5.93 


0.0000 


0.1334 


0.1332 


0.05702 


+ 11 
-12 


0.0945 


+ 4 
- 5 


5.93 


0.0000 


0.1334 


0.1334 


0.05390 


+ 12 
-12 


0.0894 


+ 4 
- 5 


5.93 


0.0000 


0.1337 


0.1337 


0.04464 


+ 12 
-11 


0.0754 


+ 4 
- 4 


5.93 


0.0000 


0.1339 


0.1337 


0.04151 


+ 12 
-11 


0.0703 


+ 4 
- 4 


5.93 


0.0000 


0.1339 


0.1339 


0.03838 


+ 13 
-12 


0.0653 


+ 4 
- 4 



TABLE XI. The values of <r(i) > and < s(t) > used to define mpcAC, see eq.(^ 



First Approach 



5.2000 0.1355 0.32 + ^ 



4 



5.2000 0.1350 0.393 +^g 

5.2600 0.1345 0.365 + « 

— 6 

5.2900 0.1340 0.349 + g 



5.9300 0.0000 0.376 +^2 



Second Approach 



0.35 t 2 



Third Approach 



0.43 t 2 



TABLE XII. J values from the various approaches as described in the text. 



o fFormil 



5.2000 0.1355 0.110 + 3 



5.2000 0.1350 0.115 '^ 3 

5.2600 0.1345 0.118 + I 

5.2900 0.1340 0.116 + I 



5.9300 Quenched 0.1186 ^H 



TABLE XIII. Lattice spacing determined from the mesonic sector using the method of 



Mps/Mv 



5.2000 0.1355 0.578 +}g 



5.2000 0.1350 0.700 +}„ 

5.2600 0.1345 0.783 + f 

— 5 

5.2900 0.1340 0.835 + ^ 



TABLE XIV. The ratio Afpg'^'^/M;;""^'^ for the dynamical data sets (i.e. with k = k^^'' = /s:™'). 



hadron 



Vector Meson 



5.2000 


0.1355 


0.449 


+ 21 

-15 


5.2000 


0.1350 


0.457 


+ 11 
-13 


5.2600 


0.1345 


0.472 


+ 7 
- 8 


5.2900 


0.1340 


0.470 


+ 15 
-15 



0.65 



0.76 
0.69 
0.66 



+ 3 

- 3 
+ 2 

- 2 
+ 3 

- 3 



5.9300 0.0000 0.475 + ^ 0.70 + 3 



Nucleon 



5.2000 


0.1355 


0.653 +J7 


1.28 


5.2000 


0.1350 


0.67 t \ 


1.32 


5.2600 


0.1345 


-X- 9 

0.72 _ 2 


1.12 


5.2900 


0.1340 


0.71 + I 


1.13 



+ 9 
-10 



+ 6 

- 5 

+ 4 

- 4 
+ 4 

- 4 



5.9300 0.0000 0.653 +^3 1.28 + \ 



Delta 



5.2000 0.1355 0.77 + \ 1.12 +" 



5.2000 


0.1350 


0.81 


5.2600 


0.1345 


0.80 


5.2900 


0.1340 


0.84 



t I 1.13 t I 

+ ^ 1.06 + ^ 



2 "■"" - 5 

2 0.95 + I 



5.9300 0.0000 0.81 + I 1.08 + t 

— 2 — 6 



TABLE XV. The fitting parameters for tlie partially quenched fit of eq.fcGh for the hadronic masses 



5.20 



0.1355 



0.13645 



+ 3 
- 3 



5.20 
5.26 
5.29 



0.1350 
0.1345 
0.1340 



0.13663 
0.13709 
0.13730 



+ 5 

- 6 
+ 3 

- 2 
+ 3 

- 3 



5.93 



quenched 



0.135202 



-11 
-11 



TABLE XVI. Values of ^crit obtained for all the datasets. 



hadron 



\ unitary 



Duiiitary 



Vector Meson 

Nucleon 

Delta 



0.476 



-18 



0.69 



0.84 



-f 2 



+ 3 
- 3 



0.66 



1.20 



0.94 



+ 9 



+ 12 



TABLE XVII. The fitting parameters for the "unitary" dataset fit of eq.(|37|) for the hadronic masses 



hadron 



An 



Ax 



Bo 



Bi 



Vector Meson 



Nucleon 



Delta 



0.492 



0.663 



+ 10 



+ 13 
-15 



0.84 



+ 2 



-0.004 



0.006 



-0.002 



+ 2 



+ 5 



0.61 



1.23 



0.91 



+ 4 



+ 6 



+ 8 
- 9 



0.015 



-0.001 



+ 1 



0.02 



+ 2 



TABLE XVIII. The fitting parameters for the combined fit of eq.(B8|) for the hadronic masses. 



1.60 
1.50 
1.40 
1.30 
1.20 
1.10 
1.00 
0.90 
0.80 







r=8a 



r=4a 



3 4 

t 



FIG. 1. Effective potential energies as a function of Euclidean time t (open symbols). The asymptotic estimates 
described in the text are shown as full symbols. The final estimated potential V{r) is indicated by the lines with 
error bands. The data correspond to (/Jjk""") = (5.20,0.1350) and r = (4a, 0, 0) (circles) and r = (8a, 0, 0) 
(diamonds). 
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FIG. 2. The static QCD potential expressed in units of ro. The dashed curve is a string model described in 
the text. 
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FIG. 3. The difference between the static QCD potential expressed in physical units and the prediction of the 
string model described in the text. For clarity, only data from the matched ensembles are shown. 
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FIG. 4. Effective mass plots for the pseudoscalar, vector, nucleon and delta for the /3 — 5.2, n^"^^ — 0.13500 
data set at k^^' = 0.13500. The horizontal lines show the fitted value for the mass (with error bars) obtained by 
the fitting approach described in the text. 
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FIG. 5. Vector niass plotted against pseudoscalar mass squared in units of ro, together with the experimental 
data points. 
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FIG. 



Vector-pseudoscalar hyperfine splitting in units of vq. 
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FIG. 7. J versus k"*^" using the approaches as described in the text. Note that the quenched data points have 
been plotted at n^'^^ = 0.132 for convenience. Approaches 2 & 3 are obtained after a chiral extrapolation and are 
shown as banded regions. The experimental value J = 0.48(2) is also shown. 



-0.05 



to 




-0.15 



-0.2 



-0.25 



2 3 

l/{aMpg(unitary)} 



FIG. 8. Si as a function of (1/Mp^) for i = y'a and the mass pairs {Mk*, Mk) & {Mp, M-n). Si is defined 
in eq. (K3) with j = ro- 
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FIG. 9. The Edinburgh plot for aU the data sets. All degenerate k™ correlators have been included. The 
phenomenological curve (from [M]) has been included as a guide to the eye. 
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FIG. 10. Effective masses for the A^^ ground state on the (/3, k) = (5.20,0.13550) ensemble. 
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FIG. 11. Vacuum expectation values for Polyakov loops at various blocking levels on the 

(/3,k) = (5.20,0.13550) ensemble. 
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FIG. 12. The measured topological susceptibility, with interpolated quenched points at the same fo. The 
radius of the dynamical plotting points is proportional to fg^. The fits, independent of the quenched points, are; 
(iii) Eqn. (P) and (iv) Eqn. (Fi). 
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FIG. 13. The measured topological susceptibility. The radius of the dynamical plotting points is proportional 
to Tq^. The fits, independent of the quenched points, are: (i) Eqn. (pl[), (ii) Eqn. (|52[), (iii) Eqn. ( p3| ) and 
(iv) Eqn. (a. 



